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I.  INTRODUCTION 


The  LABCEDE  facility  is  dedicated  towards  the  study  of  the  fluore¬ 
scence  arising  from  electron  irradiated  gases.  In  this  experiment  a  pulsed, 
collimated,  mono-energetic  beam  of  electrons  is  introduced  into  a  chamber 
containing  a  constant  pressure,  flowing  gas.  Free  radicals  and  excited 
atomic  or  molecular  states  are  created  in  collisions  between  the  electron 
beam  and  gas  target  and  these  species  can  either  fluoresce  directly  or 
undergo  further  chemical  reactions  which  ultimately  result  in  fluorescence. 

In  the  LABCEDE  facility  this  irradiation  induced  fluorescence  is  observed 
along  an  optical  axis  at  right  angles  to  the  beam  direction. 

Present  studies  are  directed  towards  the  study  of  the  infrared  radia¬ 
tion  emitted  from  electron  irradiated  gas  mixtures  of  Ar/CO^.  This  emis¬ 
sion  is  monitored  through  use  of  a  Michelson  interferometer  which  is  operated 
in  a  mode*  which  allows  observation  of  both  the  spectral  and  temporal  evolu¬ 
tion  of  the  fluorescence.  Since  the  electron  beam  is  operated  in  a  pulsed 
mode  (25  msec  period,  10%  duty  cycle),  both  the  fluorescence  excitation  and 
decay  can  be  monitored. 

The  two  dominant  sources  of  radiation  observed  in  electron  irradiated 
mixtures  of  CO^/Ar  are  the  CO^V^)  band  at  4.  3  pim  and  the  fundamental 
vibration/rotation  band  of  CO  which  spans  the  spectral  region  of  4.  5  -  6.  0 
)im.  This  latter  radiation  is  found  to  correspond  to  that  of  highly  vibra- 
tionally  excited  CO  present  in  strongly  non- Boltzmann  vibrational  distri¬ 
butions.  A  preliminary  description  of  these  observations,  along  with  a 

detailed  examination  of  the  electron  beam  gas  interaction,  has  been  pre- 
2 

sented  earlier.  A  much  larger  data  base  is  examined  in  this  report  and 
absolute  values  of  the  fluorescence  excitation  efficiencies  are  provided. 

A  detailed  examination  of  the  LABCEDE  system  fie  Id- of -view,  in¬ 
cluding  the  development  of  an  absolute  calibration,  is  presented  in  Section  II, 
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The  CCh(v  )  data  base  is  considered  in  Section  III.  Topics  discussed  in 
this  section  include  vibrational  relaxation,  band  fluorescence,  self  absorp¬ 
tion  and  transmission,  and  the  absolute  electron  excitation  efficiency  for 
the  CO_(v_)  mode.  An  examination  of  the  CO  data  base  is  provided  in  Sec- 
tion  IV.  Subjects  considered  include  the  deduction  of  the  level  dependent 
vibrational  relaxation  rate  constants  for  collisions  between  CO(v)  and  CO^, 
and  the  determination  of  the  absolute  electron  excitation  efficiency  for 
vibrationally  excited  CO  in  Ar/CO^  mixtures.  Lastly  the  summary  and 
conclusions  of  the  report  are  presented  in  Section  V. 

An  examination  of  NO  vibration/rotation  band  fluorescence  occurring 

in  electron  irradiated  air  and  N^/O^  mixtures  has  also  been  performed 

during  the  present  reporting  period.  This  work  has  been  described  in  an 
3 

earlier  report  and  is  included  here  as  an  appendix. 
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II.  CALIBRATION  STUDIES 


The  first  steps  in  developing  an  absolute  intensity  calibration  for 
the  LABCEDE  system  involve  defining  the  interferometer  field- of- view 
and  the  relative  response  function  of  the  interferometer  for  radiation  occur¬ 
ring  at  different  positions  within  this  field- of- view.  The  LABCEDE  test 
chamber  is  cylindrical  with  inner  radius  of  10  cm  and  length  of  16  cm  with 
the  electron  beam  directed  along  the  cylinder  axis.  A  circular  arsenic  tri¬ 
sulfide  window,  with  radius  of  4  cm,  is  positioned  along  one  side  of  the 
cylinder  and  centered  8  cm  downstream  of  the  beam  entrance  nozzle.  The 
optical  axis  of  the  interferometer  co-incides  with  a  line  drawn  radially  from 
the  cylinder  axis  through  the  center  of  the  window.  The  distance  between 
the  cylinder  axis  and  the  entrance  slit  of  the  interferometer  is  41  cm  and 
the  effective  internal  length  of  the  interferometer  is  17  cm.  Lastly  the 
distance  between  the  cylinder  axis  and  the  window  exterior  is  13  cm. 

The  limiting  stop  of  the  interferometer  is  the  detector  (PbS  or  PbSe) 
which  is  2mm  x  2mm.  This  detector  is  positioned  immediately  beyond  a 
short  focal  length  lens.  The  optimum  relative  position  of  the  detector  is 
set  by  requiring  that  the  detector  be  filled  by  the  central  fringe  of  the  dif¬ 
fraction  pattern  resulting  from  a  point  source  positioned  along  the  optical 
axis  at  the  center  of  the  test  tank.  Such  a  setting  allows  for  optimum  inter¬ 
ferometer  resolution  for  a  given  throughput  while  minimizing  radial  dis¬ 
tortion  and  chromatic  aberrations.  This  setting  is  achieved  by  careful 
positioning  of  the  detector  relative  to  the  lens  and  appropriate  masking  of 
the  detector.  A  complication  arises  in  that  the  central  fringes  for  off-axis 
radiation  will  fall  off  detector  center  and  thus  some  of  this  radiation  will  be 
lost.  The  net  result  is  that  the  system  response  will  be  non-uniform  across 
the  field-of-view  and  this  effect  must  be  accounted  for  in  the  development  of  an 
absolute  calibration.  (Note  this  effect  may  be  ameliorated  by  positioning 
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collimating/telescoping  optics  between  source  and  interferometer.  Such 
a  system  is  presently  under  development.  ) 

The  interferometer  field- of- view  was  mapped  for  a  source  held  at 
a  distance  of  41  cm  from  the  entrance  aperture  (corresponding  to  the  dis¬ 
tance  from  the  test  chamber  axis).  The  detailed  detector  response  within 
the  field-of-view  was  measured  using  several  different  methods: 

-  by  using  a  blackbody  point  source  and  radiometrically  measuring 
relative  response  as  a  function  of  source  position  perpendicular 
to  the  optical  axis; 

-  by  using  a  blackbody  point  source  and  interferometrically  moni¬ 
toring  relative  intensities  at  several  source  positions; 

-  by  using  a  spectrally  filtered  blackbody  point  source  and  radio¬ 
metrically  and  interferometrically  monitoring  as  above; 

-  by  using  a  planar  line  source  (He  lamp)  much  larger  than  the  field- 
of-view  and  radiometrically  observing  detector  signal  as  a  func¬ 
tion  of  circular  aperture  size  at  the  source. 

All  experiments  yielded  compatible  results;  the  field-of-view  is 
grossly  circular  about  the  optical  axis  axis  and  has  a  gaussian  decay  with 
increasing  radial  position.  The  detailed  field-of-view  as  determined  by 
the  first  method  above  is  displayed  in  Fig,  1.  In  Fig.  2,  the  fraction  of 
the  total  field-of-view  response  predicted  within  radius  r  is  plotted  using  a 
radial  distribution  as  displayed  in  Fig.  1  and  assuming  cylindrical  sym¬ 
metry.  Also  plotted  are  measurements  with  apertures  obtained  by  the  last 
method  above.  Agreement  is  seen  to  be  excellent. 

Clearly,  response  non- uniformity  precludes  use  of  the  full  extended 
field-of-view  area,  17(3.  8  cm)^,  in  developing  an  absolute  calibration.  The 
area  under  the  circular  Gaussian  response  curve  may  be  used  to  determine 
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Fig.  2  Predicted  fraction  of  total  detected  intensity  contained  within  the  normalized 
radius,  r/rmax,  using  the  radial  detector  response  function  of  Fig.  1  and 
assuming  circular  symmetry.  Prediction  is  made  from  blackbody  point  source 
data,  and  is  the  maximum  off  axis  source  displacement  which  produced 

measurable  signal.  Also  displayed  are  measurements  using  apertures  in  front 


an  average  flat  response  across  the  whole  field-of-view  or  alternately  the 
response  may  be  defined  as  unitary  for  radial  positions  less  than  an  effec¬ 
tive  radius,  and  zero  beyond  this  value.  For  the  PbSe  detector  used 

for  the  above  study  r  ^  =  1.  93  cm  at  the  distance  of  the  electron  beam  from 
the  interferometer.  Note  that  the  area  of  this  effective  field-of-view  is 
about  one-fourth  (0.  26)  that  of  the  full  geometric  field-of-view  thus  exhibiting 
the  severity  of  the  effect  of  this  non-uniform  response. 

Approximation  of  the  real  Gaussian  field-of-view  by  a  flat  response 
field-of-view  with  r^^  =  1.  93  cm  is  possible  only  for  a  constant  extended 
source  which  fills  the  full  field-of-view.  In  the  present  experiments  the 
radiation  source  is  roughly  conical  in  shape  with  the  cone  symmetry  axis 
perpendicular  to  the  interferometer  optic  axis.  The  source  radiation  will 
be  peaked  along  the  horizontal  optical  center  line  and  the  radiation  will  be 
highly  radially  unsymmetric  about  the  center  of  the  detector  field-of-view. 
Furthermore  the  radiation  source  need  not  fill  the  field-of-view. 

The  electron  beam  (and  thus  fluorescence)  spatial  distribution  may 
be  modeled  using  empirical  formulas  developed  by  Center^,  For  all  hori¬ 
zontal  positions  downstream  of  beam  origin,  x,  the  electron  beam  current 
density  in  a  vertical  plane  perpendicular  to  the  beam  axis  is  a  radially  sym¬ 
metric  Gaussian  about  that  axis.  The  radius  within  which  one  half  of  the 

beam  current  is  contained,  r  grows  with  beam  propagation  through  the 

1/& 

target  chamber: 


-28  3  -18  2 

Where  A  and  B  have  values  of  4.  0  x  10~  cm  and  2.  7  x  10"  cm  respec¬ 
tively  for  an  electron  beam  accelerated  by  30  kV  propagating  through  argon, 

17  3 

N  is  the  number  density  in  the  test  chamber,  2.  9  x  10  molecules/cm  for 
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the  case  of  interest,  and  r  is  the  beam  radius  at  the  chamber  entrance. 

o 

The  beam  current  density  at  any  point  in  space  i  (x,  y,  z)  is  then  simply: 


i(x,  y,  z) 


Ie-(r/l.  2  r1/2) 


rrr 


1/2 


(2) 


Where  I  is  the  total  current  and  the  coordinate  system  is  given  in  Fig.  3. 

The  beam  profile  presented  in  the  figure  is  appropriate  for  30  kV  and  9  torr 

of  argon.  Displacements  in  y  and  z  about  the  beam  axis  are  symmetric  and 

2  2  1/2 

use  of  r  =  (y  +  z  )  is  valid.  The  detector  field-of-view  is  taken  as  per¬ 
pendicular  to  the  beam  axis  (along  z)  and  is  centered  on  the  beam  axis  at 
x  =  8  cm,  with  the  full  field-of-view  extending  from  4.  2  to  11.  8  cm  from  the 
beam  entrance  to  the  chamber.  For  the  case  of  Fig.  3  only  89.  7%  of  the 
beam  current  between  the  planes  x  =  4.  2,  1 1.  8  cm  fall  within  this  field-of- 
view. 

The  response  of  the  interferometer  to  this  non-uniform  source  may 
be  determined  by  convoluting  the  volumetric  source  function,  Eq.  (2),  over 
the  experimentally  determined  relative  response  function  displayed  in  Fig.  1. 
This  has  been  done  for  a  source  corresponding  to  30  kV  and  9  torr  Ar  and 
it  was  found  that  the  effective  system  response  was  0.  32,  i.  e.  ,  as  if  only 
32%  of  the  source  radiation  produced  detector  signal.  This  value  is  to  be 
compared  with  the  factor  of  0.  26  deduced  for  the  case  of  a  constant  source 
filling  the  field-of-view. 

It  should  be  noted  that  the  actual  beam  shape  has  not  been  measured 
in  the  present  experiment.  The  validity  of  using  Centers  formalism  has 
been  checked  by  measuring  the  variation  in  the  observed  interferometer 
signal  with  effective  test  chamber  window  radius,  i.  e,  ,  limiting  the  field  of 
view.  The  source  for  this  test  was  infrared  radiation  emanated  by  CO_ 

M 

seeded  Ar  held  at  a  pressure  of  9  torr  and  irradiated  by  30  kV  electrons. 
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The  measured  fractional  transmitted  signal  as  a  function  of  radius  was  com¬ 
pared  to  that  predicted  by  convoluting  Eq.  (2)  over  the  measured  system 
response  function  and  this  comparison  is  shown  in  Fig.  4.  As  can  be  seen 
data  and  theory  are  in  excellent  agreement  thus  providing  an  additional  degree 
of  confidence  in  the  analysis. 

The  additional  correction  developed  here  for  source  non-uniformity 
is  approximately  a  25%  effect  and  is  very  specific  to  the  choice  of  beam/gas 
conditions.  However  the  variation  in  this  correction  is  not  severe.  The 
most  pronounced  effect  of  source  non- uniformity  corresponds  to  the  beam 
acting  as  a  line  source  along  the  cylinder  axis  and  in  this  case  it  may  be 
readily  shown  that  the  effective  relative  response  would  be  0.44.  On  the 
other  hand  as  conditions  are  adjusted  such  that  the  beam  grows  more  rapidly 
the  fluorescence  becomes  more  diffuse  and  the  relative  response  approachs 
the  constant  source  value  of  0.  26.  It  is  of  course  assumed  here  that  the 
gas  does  not  act  as  a  thick  target,  i.  e. ,  that  the  electrons  are  sufficiently 
energetic  to  traverse  the  field- of- view.  Furthermore  it  has  been  assumed 
that  the  fluorescence  occurs  co-incident  with  the  beam  current  density. 

There  will  actually  be  a  finite  time  lag  between  excitation  and  fluorescence 
and  thus  the  source  of  fluorescence  can  be  more  diffuse  than  the  electron 
source,  necessitating  a  smaller  response  correction.  For  the  character¬ 
istic  times  and  gas  densities  of  interest  in  the  present  studies  excited  state 
diffusion  distances  will  be  quite  small  and  this  latter  effect  can  be  neglected. 

It  now  remains  to  develop  an  absolute  calibration  through  use  of  a 
black  body  source.  It  is  to  remembered  that  the  fluorescence  from  test 
chamber  center  traverses  a  distance  of  13  cm  prior  to  entry  into  the  room 
air  environment.  Thus  the  black  body  was  positioned  28  cm  rather  than 
41  cm  from  the  interferometer  entrance  aperture  in  order  to  properly  repro¬ 
duce  the  effects  of  atmospheric  absorotion.  The  black  body  source  used  in 
the  calibration  allowed  for  several  different  aperture  sizes,  all  of  which  were 
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Fig.  4  Predicted  fraction  of  total  electron  beam  intensity  (PbSe  2  mm  round  detector)  con¬ 
tained  within  radius  r  using  the  experimentally  measured  radial  response  function  of 
Fig.  1  viewing  the  spatially  non-uniform  electron  beam  as  in  Fig.  3.  Also  displayed 
are  measurements  taken  with  apertures  centered  on  optical  axis  (•). 


much  smaller  than  the  area  of  the  field- of- view.  The  choice  of  blackbody 

aperture  size  is  subject  to  two  opposing  constraints.  If  the  aperture  area 

is  too  small  300  K  background  radiation  dominates  the  recorded  spectrum 

and  the  calculated  responsivity  is  much  too  small.  On  the  other  hand  if 

the  aperture  area  is  too  large,  the  detector  saturates  and  response  becomes 

nonlinear  and  too  small.  This  nonlinearity  would  be  uniform  for  the  detector 

DC  level,  however,  the  interferogram  intensity  would  be  asymmetrically 

distorted.  The  actual  experimental  data  from  a  pulsed,  low  intensity  source 

is  not  restricted  by  these  constraints.  A  blackbody  aperture  of  area  1  x 
-2  2 

10  cm  was  found  to  be  the  largest  size  allowable  for  a  temperature  of 
1173  K  and  was  used  in  the  calibration. 

After  the  calibration  interferogram  is  Fourier  transformed  a  rela¬ 
tion  between  detector  signal  and  blackbody  radiation  may  be  developed,  i.  e.  , 

SBB(X)  =  a(X)  e  Nx  (3) 

where  S  (X)  is  the  detector  signal  in  mV  induced  at  wavelength  X  by  a  black- 
BB  2 

body  source  having  intensity  N  in  W/cm  -sr-|jm  and  area  A  .  The  quantity 

X  ,  BB 

a(X)  is  the  response  function  and  e  is  the  collection  efficiency  for  point 
source  radiation  positioned  along  the  optical  axis  28  cm  from  the  interfero¬ 
meter  entrance  aperture.  This  latter  quantity  may  be  related  to  €',  the 
centerline  collection  efficiency  at  a  distance  of  41  cm  (tank  center),  by  the 
inverse  square  law,  i.  e. , 


e 


(4) 


Note  that  the  internal  interferometer  length  of  17  cm  has  been  added  to  both 
distances  in  Eq.  (4) 
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If  1^  is  the  total  volume  radiation  in  Watts/sr-(jm  emanated  within 
the  field- of- view  by  a  non-uniform  source  centered  along  the  test  chamber 
axis,  then  the  relevant  detector  signal  would  be  given  by 


Ss(»  =  a(X)  e  ^ 


(5) 


where  e  is  the  average  collection  efficiency  over  the  field-of-view  specific 

to  the  source  in  question  and  T  is  the  transmittivity  of  the  chamber  window. 

K 

The  quantity  e  maybe  related  to  the  on  axis  collection  efficiency  e  by  the 
relation 


e  =  g  S  (6) 

where g  is  the  previously  determined  factor  which  accounts  for  the  non- 
uniform  response  of  the  system,  i.  e.  ,  g  =  0.  26  for  a  constant  source  or 
0.32  for  the  e-beam  case  considered  earlier. 

.  The  final  relationship  between  source  intensity  and  detector  signal 
is  arrived  at  by  combining  Eqs.  (3)  -  (6)  with  the  result 

\  -  66  NX  ABB  SSWl/k  TX  SBB(X)1  (7) 

Expression  (7)  has  been  used  in  determining  the  absolute  intensities  pre¬ 
sented  in  the  following  sections.  Note  that  the  field-of-view,  and  the  quantity 
g,  can  vary  significantly  with  detector  change  or  repositioning.  For  the  data 
to  be  discussed  the  detector  configuration  exhibited  a  field-of-view  of  radius 
3.  2  cm.  It  was  found  that  for  this  field-of-view  and  the  standard  case  of 
30  kV  electrons  injected  into  9  torr  of  argon,  the  average  normalized  response 
was  g  s  0.  38.  This  circular  field-of-view  encompasses  77.  3%  of  the  beam 
current  density  falling  between  the  planes  at  x  a  4,  8,  11.  2  cm. 


HI.  C02  INFRARED  FLUORESCENCE 

A  number  of  fluorescence  measurements  were  performed  in  room 
temperature  mixtures  of  Ar/C02»  Typical  experimental  conditions  were 
P  =  9-11  torr,  PCO2  =  0.  1  -  0.  3  torr  with  electron  beam  parameters 
V  =  30  kV  and  I~  1  mA.  In  all  cases  fluorescence  due  to  CO_(v_)  band 

L  3  _i 

transitions  was  observed  in  the  spectral  region  of  2200  -  2400  cm  .  A 
typical  emission  history  for  this  band  is  shown  in  Fig.  5.  The  electron 
beam  is  operated  in  the  pulsed  mode,  on  2,  5  msec  and  off  22.  5  msec.  The 
history  of  Fig.  5  shows  a  continuous  increase  in  the  radiation  level  while 
the  beam  is  on  followed  by  an  exponential  decay  upon  beam  termination. 

It  is  anticipated  that  the  appropriate  excitation  mechanism  producing 
this  fluorescence  is 

ke 

e  +  COz  e  +  co2tv3>  (8) 

where  the  electrons  are  dominantly  low  energy  secondary  electrons  produced 
in  ionizing  collisions  between  the  primary  electron  beam  and  neutral  gas 
target.  The  dominant  loss  mechanisms  for  CO_(v  )  are  collisional  quenching 

C02(v3)  4  Ar-»CC>2(3v2)  +  Ar  .kg  =  1.  7  x  10-15  cm3/sec  (9) 

- 14  3 

C02(v3)  +  C02  -*C02(3v2)  +  CC>2  kjQ  =  1.  1  x  10  cm  /sec  (10) 

where  the  listed  room  temperature  rate  constants  are  from  Yardley  and 
Moore^.  C02(v3)  deactivation  by  radiative  decay. 
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(11) 


is  unimportant  inasmuch  as  the  test  chamber  gas  is  severely  optically  thick 

2 

to  band  radiation  .  From  reactions  (9)  -  (11)  it  can  be  readily  shown 
that  the  total  CO_(v_)  concentration  within  the  field- of- view  while  the  beam 

j 

is  on  is  given  by  the  relationship 


Nv, =  V  Vu  - exp(-  t/Tv 


(12) 


where  is  the  global  excitation  rate  within  the  field- of- view  for  reac¬ 
tion  (8)  and 


T  =  (kg[Ar]  +  k10[CO2J) 


-1 


(13) 


After  beam  termination  at  t  =  the  concentration  will  be  given  by 


NV  =  NV  (tb}  exp  (_  (t  '  tb)/Tv  } 
3  3  3 


(14) 


Note  that  the  observed  decay  time  of  the  fluorescence  shown  in 
Fig.  5  is  1.7  msec.  This  value  is  in  excellent  agreement  with  Tv  as 
evaluted  from  Eq.  (13)  using  the  measured  rate  constants  for  k^  and 
Given  the  total  measured  band  intensity  Iv^  and  the  observed  value  of  Tv^ 
the  total  rate  of  electron  excitation  Rv^  may  be  deduced  from  Eq.  (12) 
through  use  of  the  relationship 


I  =  N  hv  £/(4tt  t  _  ) 

v3  v3  R' 


(15) 
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where  hv  is  the  energy  of  a  band  photon  and  f  is  the  specific  transmit¬ 
tance  of  the  COJv,)  band  along  the  optical  path. 

The  ultimate  goal  of  this  analysis  is  to  determine  the  absolute  ef¬ 
ficiency  of  the  electron  beam  for  promoting  CO_(v,,)  excitation.  As  pre- 
2  ci 

viously  discussed  this  may  be  deduced  directly  from  the  quantity  R-^. 

The  prime  difficulty  in  this  analysis  is  in  evaluating  the  quantity  f,  inasmuch 
as  the  COJv.)  fluorescence  is  not  only  severely  self  absorbed  within  the 

C  J 

test  chamber,  but  is  also  affected  by  absorption  by  atmospheric  CO^  along 
the  45  cm  optical  path  from  the  chamber  window  to  the  detector.  The  major, 
portion  of  this  section  is  dedicated  to  a  careful  evaluation  of  the  quantity  f. 

A.  The  Transmission  Function  for  an  Individual  Transition 

Previously  ,  the  transmitted  emission  from  the  CO_(v_)  band  was 

“  J 

calculated  considering  the  spectral  lines  to  be  exclusively  Doppler  or 
Lorentz-collisionally  broadened.  It  was  noted,  however,  that  the  actual 
shape  of  the  fluorescent  lines  has  contributions  from  both  broadening  mech¬ 
anisms,  and  that  the  transmission  losses  in  the  chamber  for  the  real  line- 
shape  would  be  intermediate  to  the  two  previous  calculations.  Accordingly, 
in  the  present  work,  a  transmission  algorithm  has  been  developed  for  the 
Voigt  profile  and  used  to  predict  the  spectral  emission  at  the  chamber 
window.  Furthermore  in  the  earlier  analysis  the  effect  of  atmospheric  CO^ 
absorption  of  the  fluorescence  was  taken  into  account  by  means  of  a  low- 
resolution  correction  to  the  experimental  data,  which  will  be  described  more 
fully  later.  In  the  present  work,  this  correction  was  not  applied  to  the  data, 
but  rather,  the  predictions  were  corrected  to  include  this  effect.  As  a  result, 
the  CC>2  spectral  data  presented  here  will  be  different  in  appearance  than  the 
"corrected"  experimental  data  presented  previously.  The  uncorrected  data 
in  both  cases  is  very  similar  in  appearance.  Thus,  the  present  approach 


is  to  present  the  experimental  spectrum  obtained  from  the  on-detector 
intensities  and  compare  it  with  predictions  which  are  corrected  for  the 
various  transmission  losses. 

From  the  electron  beam  shape  shown  in  Fig.  3,  it  is  seen  that  the 
radiating  gas  will  be  mainly  localized  in  the  field- of- view.  For  the  present 
predictions,  the  simplifying  assumption  has  been  made  that  the  vibration- 
ally  excited  gas  is  located  in  a  plane  perpendicular  to  the  detection  system 
optical  axis  at  the  center  of  the  chamber.  The  fluorescence  must  then 
travel  on  average  through  11  cm  of  room  temperature  vibrationally  equili¬ 
brated  gas  prior  to  reaching  the  chamber  window. 

The  radiation  per  unit  length  from  an  optically  thin  source  with  an 
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emitter  concentration  of  NR  molecules  cm  is 

Nr  Ibb  W/cm3-sr-cm"1  (16) 

where  IRB  is  the  blackbody  radiation  function  and  k  is  the  absorption  coef¬ 
ficient  at  frequency  V  for  the  appropriate  lineshape.  Attenuation  of  this 
radiation  over  a  distance  L  is  given  by 

i  -  \  nr  *bb  ex*(-  KNL)  (17) 

where  K  is  the  total  absorption  coefficient  at  frequency  v  due  to  all  transi¬ 
tions  (of  all  absorbers  present)  which  overlap  the  emission,  i.  e. , 

KN  =  E  k  N.  (18) 

i  V.  1 

The  fractional  transmission  of  a  single  fluorescence  line  is  given  by 
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Ibb  u(v)dv 


(19) 


where  u(v)  is  the  functional  form  of  the  absorption  coefficient  frequency 
dependence  relative  to  its  centerline  value. 


Evaluation  of  the  transmission  function  for  the  band  requires  know¬ 
ledge  of  all  ky.  at  frequency  intervals  small  compared  to  the  molecular  line- 
width  and  the  solution  of  Eq.  (19)  at  a  large  number  of  frequencies  for  each 
emission  line.  This  rigorous  approach  is  both  costly  and  time  consuming. 
However,  if  the  absorption  of  each  emitted  line  is  assumed  to  occur  only 
by  the  same  transition  which  gave  rise  to  the  radiation,  then  the  sum¬ 
mation  of  Eq.  (18)  is  not  required,  considerably  simplifying  the  calculation. 
In  this  limit  because  the  blackbody  function  is  essentially  constant  over  the 
molecular  linewidth,  Eq.  (19)  reduces  to 


T 


CD 


u(v)  exp  (-  kvNR 


L)dv 


u(v)dv 


(20) 
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In  the  list  report,  this  transmission  function  was  evaluated  for 
CO^(v^)  band  transitions  for  both  Doppler  and  Lorentz  broadened  lineshapes. 
The  enhanced  emission  in  the  wings  of  the  collision-broadened  lineshape 
resulted  in  a  significantly  different  fractional  transmission  than  for  Doppler 
lines,  about  a  factor  of  five  greater  for  a  chamber  pressure  of  1  torr  CO^, 
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9  torr  argon.  Because  of  this  large  difference,  development  of  a  more 
detailed  transmission  prediction  for  a  lineshape  more  closely  approxi¬ 
mating  the  true  emission  lineshape  was  warranted. 

Doppler  and  collisional  broadening  are  entirely  independent  pro¬ 
cesses  and  a  composite  lineshape  in  which  both  contribute  may  be  viewed 
as  one  in  which  each  frequency  element  of  the  inhomo geneou sly  broadened 
collisional  (Lorentz)  distribution  is  broadened  by  a  Doppler  distribution. 
The  resultant  lineshape,  referred  to  as  a  Voigt  profile,  is  defined  by 


k 

V 


-«gb_r2L. 
2  2 
a  +(S-y) 


dy 


(21) 


where  kv  is  the  centerline  absorption  coefficient  for  a  Doppler  broadened 
CL 

line  with  a  linewidth  of  b^  (half  width  at  half  maximum,  HWHM)  which  may 
be  specified  by  the  relation 


k  =  -5-  jLni/n'  (22) 

VC,  bD 

L 

where  S  is  the  integrated  line  strength  for  a  single  vibration  rotation  transi¬ 
tion  in  units  of  cm/molecule.  The  quantity  a  in  Eq.  (6)  is  the  ratio  of  the 
Doppler  and  collisional  half  widths 


(23) 


y  is  specified  by 
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(24) 


(V  -vc  ) 


and  §,  which  represents  the  doppler  lineshape  variation  with  frequency,  is 
defined  as 


S=^(V-VC  )  (25) 

bD  CL 

t 

Integration  over  the  dummy  frequency  V  ,  permits  the  inclusion  of  contribu¬ 
tions  from  all  Doppler- shifted  transitions  in  the  total  absorption  coefficient 
at  frequency  v. 

The  room  temperature  Doppler  broadening  of  the  CO.  molecular 
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lines  at  the  transition  frequency  is  2.  2  x  10  cm  HWHM.  Collisional 

broadening  of  the  lines  will  include  contributions  from  Ar  at  9  torr  (5.  8  x 

10  ^  cm  *  HWHM)  and  CO^  self-broadening  (1.  3  x  10-5  -  1.  3  x  10~^  cm’' 

HWHM)  for  0.  1  -  1.0  torr  of  CC^.  Their  sum  is  the  collisional  half-width, 

b  ,  for  a  given  experiment, 
c 

Equation  (21)  provides  the  functional  form  of  u(v)  and  substitution 
into  Eq.  (20)  yields 


which  is  the  fractional  transmission  for  an  optically  thin  Voigt-broadened 
emission  source  attenuated  by  resonant  non- overlapping  Voigt  absorbers. 
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Direct  evaluation  of  Eq.  (26)  would  require  integration  over  the 

i 

dummy  variable  v  at  many  values  of  frequency  v  over  each  rota¬ 
tional  line.  Numerical  tabulations  of  exact  Voigt  profiles  exist  but  closed- 
form  approximations  which  permit  calculations  of  the  line  shape  are  more 


generally  applicable.  Several  algorithms  have  been  developed  to  permit 

7-9  7 

numerical  evaluation  of  for  Voigt  lines.  Whiting’s  approximation  to 
the  true  line  shape  has  been  applied  to  analyze  the  present  experiments. 


The  Voigt  linecenter  intensity  was  first  calculated  as 


k  - - § - T  (27) 

CT  2b  [1.065  +  0.447d  +  0.058  d  ] 

L  v 


where  d  is  the  ratio  of  the  Voigt  and  collisional  halfwidths,  b^/b^,  and 


'.v=|[bc  +  (b%4^,1/2]  . 


Equation  (27)  is  analogous  to  Eq.  (22)  which  is  for  Doppler  lines.  Equations 
(27)  and  (28)  have  been  found  to  be  accurate  to  within  0.  5%  of  the  exact  Voigt 
linewidths  and  centerline  intensity.  The  Voigt  absorption  coefficient  at  any 
frequency  v  relative  to  line  center  is  approximated  by  the  empirical  relation 
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where  z  is  (v-  v<3  ).  This  expression  has  been  determined  to  be  accurate 

J-J 

to  within  2.  5%  of  the  exact  Voigt  profile  for  all  v  and  is  an  empirical  repre¬ 
sentation  of  Eq.  (21).  Substitution  of  Eqs.  (27)  -  (29)  into  Eq.  (26)  results 
in  a  considerable  simplification,  i.  e. , 


T 


CD 


k  /S  exp[-(SN  L)(k  /S)]dy 

V  K  V 


CD 

I 

-00 


kv/S  dy 


(30) 


The  fractional  transmission  of  a  single  rotational  line  over  its  entire 
bandpass  may  now  be  calculated  as  a  function  of  S,  the  linestrength.  A  plot 
of  t  vs.  SNj^L  using  Whiting's  approximation  for  the  true  Voigt  line  shape  is 
given  in  Fig.  6.  In  addition  to  the  Voigt  fractional  transmission,  the  pre¬ 
viously  derived  Doppler  and  Lorentz  transmission  functions  are  also  plotted. 

It  is  seen  that  the  Voigt  transmission  shares  an  asymptote  with  the  Doppler 
prediction  for  small  values  of  SN_  L.  At  large  values  of  SN  L,  the  wings 

£v  x\. 

of  the  lineshape  become  important  and  the  Lorentzian  component  of  the  Voigt 
profile  dominates  the  dependence  of  the  fractional  transmission  on  line- 
strength.  These  transmission  predictions  were  used  in  conjunction  with 
molecular  linestrength  calculations,  which  will  be  described  next,  to  generate 
predicted  emission  spectra  at  the  chamber  window. 

B.  Transmission  of  the  CO^v^)  Band  Within  the  Test  Chamber 

The  frequencies  and  linestrengths  were  calculated  for  rotational 
transitions  of  the  fundamental  asymmetric  stretch  of  the 
(626)  molecules.  Also  included  in  the  synthesized  spectrum  were  hot-band 
transitions  with  Av3  =  1  involving  all  significantly  thermally  populated 
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Fractional  transmission  as  a  function  of  rotational  linestrength  times  column 
density  for  Doppler,  Voigt,  and  Lorentz  calculations.  P.  =  9  torr 


vibrational  modes,  as  well  as  the  fundamental  and  hot-band  transitions  of 

,  -.13  16  _  16  12  18  .  16  _  12  17  .  .  .  .  . 

the  C  O  ,  O  C  O  ,  and  O  C  O  isotopic  species.  The  band 

center  frequencies  and  integrated  vibrational  band  intensities  are  given  in 

Table  I  for  all  the  transitions  included  in  the  calculation.  The  detailed 

spectroscopic  modeling  of  the  rotational  line  spacings  and  degeneracies 

2 

involved  in  the  spectral  synthesis  has  been  described  previously  and  will 
only  be  outlined  here.  The  constants  on  which  the  calculations  are  based 
were  taken  from  the  AFCRL  Atmospheric  Line  Parameters  Compilation**. 

The  total  integrated  bandstrength,  S^,  is  partitioned  among  the 
allowed  rotational  transitions  according  to  the  formula 


S 


v  J 

—  exp- 
rot 


(31) 


where  S  is  the  individual  rotational  linestrength,  Q  .is  the  rotational 

f 

partition  function,  J  is  the  rotational  level  from  which  the  transition 

originates,  and  S  are  the  linestrength  factors  given  by  the  Honl- London 
12  ^ 

formulae.  The  emitted  intensity  per  molecule,  I,  can  then  be  defined  for 
each  rotational  line  by  integrating  Eq.  (16): 


I  = 


00 

J  <VNR>  dv  * 

JCO 


dv  =  I  S  (W/molecule-sr) 


(32) 


where  the  relation  5=  /  k  d  has  been  utilized.  The  fractional  trans- 

■00  ^  ^ 

mission,  t  of  each  line  through  the  Voigt  broadened  absorbing  gas  was 
calculated  from  Eq.  (30)  over  the  range  of  0.  0  -  1.0  torr  CO^  pressure. 
The  total  emitted  intensity  at  the  chamber  window  was  obtained  as  the  sum 
of  contributions  from  each  rotational  line  from  all  the  bands. 
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TABLE  I 


CC>2(V3)  SPECTRAL  PROPERTIES 


Isotope 

— 

Transition 

Ground 

State 

V1V2  V3 

* 

s 

V 

cm/molec. 

vband* 

-1 

cm 

Branches 

Rotational  Levels 

J 

626 

00°0 

9.  598E-17 

2349,  146 

P,R 

even  J  only 

626 

o 

o 

H 

o 

7.  66E-18 

2336.  637 

P,  Q.R 

all  J,  c  odd 

636 

00°0 

9.  602E-19 

2283. 490 

a 

P.R 

even  J 

628 

00°0 

3.  33E-19 

2332.  112 

P,R 

all  J 

626 

02Z0 

3.  08E-19 

2324. 148 

P,Q,R 

all  J,  c  even 

626 

02°0 

1934E-19 

2327. 432 

P,R 

even  J 

626 

10°0 

1.  183E-19 

.  2326.  594 

P.R 

even  J 

636 

01*0 

8.  18E-20 

2271. 763 

P,Q,R 

all  J,  c  odd 

627 

00°0 

6.  37E-20 

2340.01 

P.R 

all  J 

628 

O 

H 

H 

o 

2.  58E-20 

2319. 738 

P#  Q»R 

all  J,  c&d  each 

626 

0330 

1.22E-20 

2311. 675 

P,  Q.R 

all  J,  c  odd 

627 

O 

H 

H 

o 

4.97E-21 

2327.  575 

P,  Q,  R 

all  J,  c&d  each 

636 

02Z0 

3.  46E-21 

2260.  045 

P,  Q,R 

all  J,  c  even 

This  prediction  for  the  total  transmitted  intensity  per  molecule  is 

plotted  in  Fig.  7  as  a  function  of  CO.  p  -essure.  Also  displayed  are  the 

^  2 

predictions  from  previous  calculations  for  both  Lorentz  and  Doppler 
broadened  absorbers.  As  expected  the  predicted  Voigt  intensity  falls 
between  the  values  of  the  two  earlier  calculations  at  low  CO^  pressures, 
but  follows  more  closely  the  pressure  dependence  of  the  Lorentzian  case. 

The  lower  solid  curve  in  the  figure  will  be  discussed  shortly.  As  required 
in  the  limit  of  no  attenuation  (very  low  CO^  pressure), 

E  E  I  T  =  E  S  I  (v  )  (33) 

j  i  1J  j  vj  BB  bj 

where  the  summation  over  i  includes  all  rotational  lines  and  j  includes  all 
transitions  considered. 

The  relative  contributions  of  the  hot  bands  and  isotopic  transitions 

included  in  the  calculation  is  demonstrated  in  Figs.  8  and  9.  In  Fig.  8,  the 

fractional  Voigt  transmissions  (E  T..)  are  shown  as  a  function  of  CO  pres- 

i  ij  2 

sure  for  four  of  the  strongest  transitions  in  Table  I:  the  626  fundamental 

1  2 

and  first  and  second  excited  bending  mode  transitions  (01  1  and  02  1)  and 
the  636  isotopic  fundamental.  Although  all  transitions  suffer  attenuation 
which  increases  with  CC>2  pressure,  the  626  \>3  fundamental  band  radiation 
will  be  very  severely  attenuated,  only  3%  transmitted  at  1  torr  pressure. 

The  weaker  bands  (not  shown)  will  undergo  only  slight  self-absorption  and 
thus  will  be  enhanced  relative  to  the  fundamental. 

Accordingly,  even  though  the  626  fundamental  contains  93%  of  the 
total  bandstrength  of  all  the  transitions  (as  given  by  the  right  hand  side  of 
Eq.  (33)),  it  will  be  responsible  for  only  about  50%  of  the  emission  inci¬ 
dent  on  the  chamber  window  as  is  seen  in  Fig.  9.  The  predicted  fraction  of 
the  total  radiation  intensity  which  is  due  to  each  band  is  plotted  in  this  figure 
as  a  function  of  COz  pressure  for  the  four  bands  listed.  The  strongest  hot 
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TOTAL,  TRANSMITTED  INTENSITY  (W/MOLECULE  -  STR) 


LORE 


mii 


C02  PRESSURE  (torr) 


Fig.  7.  Total  transmitted  intensity  for  all  bands  as  a  function  of  CO 
pressure  for  Yoigt  (solid),  Doppler  and  Lorentz  (broken) 
transmission  functions.  The  lower  solid  curve  is  the  Yoigt 
prediction  corrected  for  atmospheric  CO-  attenuation. 

PAr  =  9  torr’  L  =  11  cm- 


COz  PRESSURE  (torr) 


Fig.  8 


.  Fractional  transmission  of  four  of  the  strongest  bands 
of  COz  pressure  for  the  Voigt  -  broadened  absorbers. 
L  ■  1 1  cm. 


as  function 
PAr  =  9  torr. 


Fig.  9.  Relative  contributions  of  four  of  the  strongest  bands  to  the  total 
transmitted  intensity  at  the  chamber  window  as  a  function  of  CO, 
pressure.  P  =  9  torr.  L  =11  cm. 
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band  transition  is  greatly  enhanced  relative  to  its  unattenuated  value,  com¬ 
prising  one-quarter  to  one- third  of  the  total  radiation  over  most  of  the 
pressure  range.  It  is  seen  that  even  hot-band  self-absorption  becomes 
significant  at  pressures  greater  than  0.  2  torr, 

C.  The  Effect  of  Atmospheric  CO^ 

The  absorption  by  ambient  CO^  present  in  the  laboratory  air  in  the 

45  cm  optical  path  from  the  window  to  the  detector  was  taken  into  account 
2 

previously  by  simply  modifying  the  experimental  spectral  data  using  a 

blackbody  spectrum  which  included  this  atmospheric  absorption  feature. 

This  correction  (with  10  cm  *  resolution)  operated  more  or  less  uniformly 

on  all  the  emitted  lines.  In  reality,  however,  absorption  of  the  fluorescence 

exiting  the  chamber  window  by  Lorentz-broadened  room  air  CO^,  will  occur 

only  when  spectral  overlap  occurs.  The  collisional  broadening  of  CO^  by 

nitrogen  has  been  measured*^  as  0.  07  cm"*  atm-^,  Only  fluorescence 

lines  occurring  within  approximately  an  cm'*  of  strong  atmospheric  lines  are 
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significantly  absorbed.  Because  the  Voigt  linewidths  (b^  ~  2.  5  x  10  cm 
HWHM)  of  the  CO^  fluorescence  are  much  narrower  than  the  atmospheric 
lines,  calculation  of  the  atmospheric  absorption  of  *his  emission  is  greatly 
simplified  since  the  Lorentzian  line  shape  of  the  absorber  does  not  vary 
significantly  over  the  width  of  the  Voigt  line.  Indeed  at  the  off-line  center 
frequency  corresponding  to  the  maximum  slope  of  the  Lorentz  lineshape, 
for  the  Lorentz  line  changes  by  about  1%  over  the  frequency  range  cor¬ 
responding  to  8  b^  (1%  intensity  points  of  the  Voigt  profile). 

Consequently,  the  emission  lines  may  be  considered  as  delta  func¬ 
tions  in  frequency  space,  and  the  atmospheric  absorption  determined  as  the 
amount  of  attenuation  of  each  emission  line  by  all  absorbing  atmospheric 
transitions.  The  effects  of  overlapping  lines  and  non- resonant  absorption 
maybe  included  in  this  rigorous  (Hi-tran)  calculation.  Although  the  dominant 
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atmospheric  absorbing  transition  is  the  626  fundamental,  the  first  hot  band 
and  636  fundamental  transitions  make  significant  contributions  to  the  absorp¬ 
tion.  Inclusion  of  the  remaining  10  atmospheric  CO^  absorption  bands  in 
Table  I  produced  a  negligible  effect  on  the  predictions.  The  appropriate 
rotational  linestrengths  and  shapes  were  used  in  the  calculation  for  each 
CC>2  transition.  The  Lorentz  broadened  lineshape  is  given  by 


k 

V 


o  (v-v  ) 

o 


where  k  is  the  line  center  absorption  coefficient 
o 


(34) 


(35)  . 


If  the  concentration  of  atmospheric  CO_  is  taken  as  the  standard  mixing  ratio 

2  17  -2 

(0.  033%),  the  resultant  optical  path  column  density  is  3.  5  x  10  cm  . 

Atmospheric  attenuation  of  the  fluorescence  was  calculated  for  this  column 
density  over  a  range  of  CO^  chamber  pressures.  The  total  transmitted 
intensity  (W/molecule- str)  for  all  rotational  lines  of  all  the  CO^  vibra¬ 
tional  bands  is  plotted  as  the  lower  solid  line  in  Fig.  7.  This  curve  includes 
the  effects  on  the  emission  of  self- absorption  by  Voigt  broadened  lines  in  the 
chamber  (upper  solid  line)  and  the  subsequent  atmospheric  COg  attenuation, 
and  thus  represents  the  prediction  of  radiation  per  molecule  which  will  reach 
the  detector.  The  atmospheric  attenuation  is  a  relatively  constant  factor  of 
two  over  most  of  the  chamber  CO ^  pressure  range  considered.  Furthermore 
the  predicted  intensity  per  molecule  reaching  the  detector  varies  by  less 
than  a  factor  of  five  over  the  chamber  CO£  pressure  range  0-1  torr. 


-33- 


Relative  contributions  of  the  three  strongest  bands  to  the  total  pre¬ 
dicted  intensity  at  the  detector  is  plotted  in  Fig.  10  over  the  range  0.  0  - 
0.  5  torr  CO_  pressure.  Tne  626  fundamental  radiation  is  severely  attenuated 
over  the  45  cm  optical  path  in  the  room;  so  much  so,  that  the  hot  band  (which 
has  a  band  strength  which  is  only  8%  of  the  fundamental’s)  is  the  dominant 
contributor  to  the  total  radiation  intensity  at  nearly  all  CO^  pressures  con¬ 
sidered.  Comparison  of  Figs.  9  and  10  illustrates  this  effect  clearly.  Even 
the  636  isotope  fundamental  transition  will  be  responsible  for  more  of  the 
total  transmitted  radiation  than  the  626  fundamental  at  pressures  above 
0.  3  torr. 

Figure  10  is  basically  similar  to  the  analogous  Doppler  predictions 

2 

performed  previously  which  did  not  take  into  account  atmospheric  C02 
attenuation.  It  may  be  seen  from  Fig.  7  that  the  previous  predicted  total 
intensity  for  Doppler  lines  differs  by  less  than  50%  from  the  present  more 
rigorous  calculation  over  all  but  the  lowest  pressures  considered.  Figure  7 
provides  proof  of  the  previous  conjecture  that  the  effects  of  Voigt  absorbers 
and  ambient  CO^  would  be  offsetting,  and  that  Doppler  profile  predictions 
provided  a  reasonable  estimate  of  the  total  transmitted  radiation. 

The  above  transmission  predictions  are  based  on  the  CO^  concentra¬ 
tion  in  the  optical  path  being  equal  to  the  standard  atmospheric  mixing  ratio. 
In  an  attempt  to  minimize  the  effects  of  atmospheric  CO^,  efforts  were  made 
to  purge  the  optical  path  of  CO^.  The  degree  of  success  will  be  evaluated 
later  from  spectral  shape  predictions  and  blackbody  absorption  features. 

The  effect  of  changes  in  the  atmospheric  CO^  concentration  on  the 
total  transmitted  intensity  predictions  in  shown  in  Fig.  11.  The  fraction  of 
chamber  window  fluorescence  transmitted  to  the  detector  is  plotted  as  a 
function  of  chamber  CO^  pressure  for  a  number  of  atmospheric  CO^  concen¬ 
trations  over  the  range  33  -  660  ppm.  (The  standard  atmospheric  mixing 
ratio  is  330  ppm. )  Note  that  at  the  lower  test  chamber  CC>2  pressures,  the 
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FRACTION  OF  THE  TOTAL  TRANSMITTED  INTENSITY  AT  THE  DETECTOR 


l 


EXPERIMENTAL  C02  PRESSURE  (torr) 

Fig.  10.  Relative  contributions  of  the  three  strongest  bands  to  the  total 
transmitted  intensity  after  attenuation  by  both  chamber  and 
atmospheric  C02,  plotted  as  a  function  of  chamber  C02  pressure. 
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FRACTION  TRANSMITTED  BY  ATMOSPHERE 


EXPEIUMENTAL  COz  PRESSURE  (torr) 

fig.  II,  The  fraction  of  the  total  COg  fluorescence  at  the  chamber  window 
transmitted  through  an  atmospheric  pressure  optical  path  con¬ 
taining  various  concentrations  of  CO,  plotted  as  a  function  of  the 
experimental  COg  pressure, 

-36- 


chamber  window  fluorescence  is  dominated  by  the  626  fundamental  radia¬ 
tion,  which  is  severely  resonantly  absorbed  by  the  atmospheric  CO^.  As 
a  result,  the  fractional  transmission  of  the  chamber  window  fluorescence 
is  lowest  for  this  extreme.  Nonetheless  for  CO^  pressures  >  0.  1  torr  the 
fractional  transmission  varies  by  only  a  factor  of  2  -  2.  5  over  the  full  range 
of  atmospheric  CO^  concentrations  shown  and  also  the  total  additional 
attenuation  is  less  than  a  factor  of  three.  The  total  intensity  incident  on 
the  detector  (W /molecule -str)  for  a  given  atmospheric  CO^  pressure  is 
calculated  as  the  product  of  the  attenuation  of  Fig.  11  and  the  Voigt  chamber 
window  predictions  of  Fig.  7. 


D.  Comparison  With  the  Data  Base 

In  order  to  permit  comparison  with  the  experimental  spectra,  the 
individual  transmitted  line  intensities  must  be  convolved  with  a  system  slit 

t 

function,  g(v,  V  ).  The  appropriate  function  for  the  LABCEDE  interfero¬ 
metric  detection  system  is  a  sine  function  defined  by 


g(v,v') 


0.  82  Av 

i 

IT  (V-V  ) 


(36) 


which  determines  the  contribution  at  frequency  v  from  a  line  at  frequency 

i 

v  for  a  system  resolution  Av.  For  all  experimental  runs  sufficient  inter- 
ferogram  points  were  taken  to  yield  10  cm  *  resolution  spectra.  Because 
the  fluorescence  lines  are  much  narrower  than  the  system  resolution,  only 
their  line  center  separation  from  the  convolution  frequency  is  required. 

The  sine  function  is  truncated  at  the  second  zero  crossing  on  each  side  of 
the  center  frequency.  All  the  transmission  predictions  for  the  individual 
V-  band  transitions  (either  at  the  chamber  window  or  at  the  detector)  were 

*  _  i 

sorted  by  frequency,  summed  in  0.  4  cm  cells  and  convolved  with  the  sine 
function  to  yield  a  predicted  spectrum  for  comparison  with  the  experimental 
data. 
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Spectra  from  sixteen  runs  were  included  in  the  analysis.  The  rele¬ 
vant  experimental  parameters  for  these  runs  are  listed  in  Table  IL  The 
CC>2  pressures  of  these  runs  were  kept  in  the  0.  1  -  0.  2  torr  range  in  order 
to  permit  observation  of  CO  relaxation  kinetics  (see  Section  IV).  The 
observed  CO^  fluorescence  near  beam  termination  from  a  slowly  purged 
runs  is  shown  in  Fig.  12.  The  dominant  spectral  feature  is  the  broad  peak 
at  *-2310  cm  *  with  a  shoulder  at  slightly  higher  frequency.  A  second  strong 
feature  at  2350  cm”*  and  the  smaller  peak  at  ~  2380  cm"*  were  also  very 
reproducible  from  run  to  run.  The  shape  of  the  broad  feature  centered  at 
2270  cm"*  varies  from  run  to  run,  but  its  relative  intensity  remained 
constant. 

The  CC^  fluorescence  spectra  for  four  runs  with  very  similar  experi¬ 
mental  and  purge  conditions  (11MAR2,  16MAR2,  16MAR3,  and  16MAR4) 
are  all  plotted  in  Fig.  13  normalized  to  the  peak  of  the  2310  cm  feature  to  . 
demonstrate  data  reproducibility.  For  comparison,  the  chamber  window 
fluorescence  predictions  for  Voigt  broadened  lines  are  shown  in  Fig.  14  -  17 
for  experimental  CC>2  pressures  of  0.  33,  0.  2,  0.  10  torr  and  trace  amount 
of  CO2  respectively.  It  can  be  seen  by  comparison  with  Fig.  12  that  the 
Voigt  predictions  without  atmospheric  attenuation  are  very  different  from 
the  data.  Figure  17  corresponds  to  the  prediction  for  the  case  of  an  optically 
thin  gas,  i.  e. ,  no  self  absorption.  In  this  case,  the  spectrum  is  primarily 
due  to  the  P  and  R  branches  of  the  626  fundamental  transition.  As  the 
absorber  concentration  increases,  other  transitions  conbribute  significantly 
to  the  spectrum  as  shown  in  Fig.  9.  The  broken  line  curve  in  Fig.  14 
represents  the  predicted  spectrum  when  only  the  626  fundamental  transition 
is  included  in  the  calculation,  clearly  demonstrating  that  the  hot  band  and 
isotopic  radiation  occurs  at  lower  wavenumbers  (to  the  red  of  the  fundamental). 

The  corresponding  spectral  predictions,  which  include  the  effect  of 
atmospheric  CC>2  attenuation,  are  displayed  in  Fig.  18  -  21.  For  these 


TABLE  II 


RUNS  INCLUDED  IN  COz  SPECTRAL  ANALYSIS 


NORMALIZED  INTENSITY 


1.  00 


Fig.  12  Normalized  experimental  data  from  11  MARCH  RUN  Z, 
0.  1  msec  before  beam  termination.  =  1  torr, 

?Ar  =  9.  1  torr,  I  =  0.81  mA,  V  =  30  kV.  Maximum 
intensity  =  3.  7  x  10"6  W/str-cm-1. 
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1.  00 


WAVENUMBER  (cm"1) 


Fig.  13  Data  from  11  MARCH  RUN  2,  16  MARCH  RUN  2,  16  MARCH 
RUN  3,  and  16  MARCH  RUN  4,  0.  1  msec  before  beam  termi¬ 
nation,  normalized  to  the  same  scale.  P^q-  =  *  torr, 

PAr  =  9  torr,  1*0.81  +  0.03  mA,  V=30  kV. 
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Fig.  14  Predicted  fluorescence  spectrum  transmitted  to  the  chamber 
window  through  Voigt  broadened  absorption  lines.  Resolu¬ 
tion  is  10  cm'l  Pqq^  =  ^  Broken  curve  is  pre¬ 

dicted  fluorescence  when  only  the  626  fundamental  is  included 
in  the  calculation.  Maximum  intensity  =  2.  28  x  10" 26  y \f  / 
molecule  -  st  r- cm" 
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Fig.  15  Predicted  fluorescence  at  chamber  windows  for  Voigt  lines 
with  10  cm'  *  resolution.  PcO£  =  0*20  torr.  Maximum 
intensity  =  2.91  x  10“  26  W /molecule-  str- cm' 
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Fig. 


17  PrIdl*ed  fI?orescence  at  chamber  window  for  Voigt  lines 
with  10  cm-1  resolution.  Trace  amount  of  CO  5L  ” 
intensity  =  4.  39  x  NT*®  W/molecuHrl  if 
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Fig.  18  Predicted  fluorescence  spectrum  at  the  detector  after  trans 
mission  through  chamber  and  atmospheric  absorbers. 
Resolution  =  10  cm' ^ =  33  torr.  Maximum  inten¬ 

sity  =  1.  25  x  10'^k  w/molecule-str-cm"^. 


NORMALIZED  INTENSITY 


Fig.  19  Predicted  spectrum  at  detector  including  effects  of  both 
chamber  and  atmospheric  absorbers.  Resolution  = 

10  cm"*.  ^C02  =  0*  2  torr.  Maximum  intensity  = 

1.  54  x  10"^  W/molecule-str-cm"*. 
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predictions  the  atmospheric  mixing  ratio  of  0.  033%  was  assumed.  All 
spectral  shapes  are  greatly  modified  from  those  shown  in  Figs.  14  -  17 
even  though  the  total  transmitted  intensity  (W/molecule- str)  changes  by 
less  than  a  factor  of  3  for  the  0.  1,  0.  2,  and  0.  33  torr  cases.  Note  that  a 
large  fraction  of  the  transmitted  intensity  occurs  in  the  hot  band  and 
isotope  emission  region  as  predicted  in  Fig.  11.  The  predictions  of 
Figs.  19  and  20  are  very  similar  in  appearance  to  the  data  of  Fig.  12, 
which  has  been  reproduced  in  Fig.  20  as  a  broken  line  curve.  The  fre¬ 
quencies  and  shapes  of  all  spectral  features  are  matched  very  well, 
although  the  relative  intensity  of  these  features  are  not  predicted  as  well. 

Because  the  optical  path  was  being  purged,  the  amount  of  atmos¬ 
pheric  CO^  present  was  uncertain.  Consequently,  the  amount  of  CC>2 
included  in  the  predictions  was  varied  in  an  attempt  to  better  match  the 
spectral  data.  Predicted  spectra  for  the  case  of  0.  1  torr  CC>2  chamber 
pressure  and  a  range  of  atmospheric  COz  concentrations  are  displayed  in 
Fig.  22.  Both  the  relative  spectral  intensity  and  the  shape  of  the  predicted 
radiation  change  over  this  range.  Most  notably  the  ~  2310  cm"  *  feature 
increases  in  relative  intensity,  narrows,  and  shifts  to  the  red  as  the  ambient 
CC>2  concentration  is  increased.  The  best  spectral  fit  was  determined  to  occur 
when  the  atmospheric  CO^  concentration  was  taken  to  be  0.  9  that  of  the 
standard  mixing  ratio,  300  ppm.  A  comparison  of  the  prediction  for  this 
mixing  ratio  and -the  experimental  spectra  of  Fig.  12  are  shown  in  Fig.  23. 

The  position  of  the  peak  and  shoulder  of  the  2310  cm’*  feature  are  well 
matched  but  the  intensity  of  the  shoulder  is  underestimated.  The  two  peaks 
at  high  wavenumber  s  are  still  overestimated.  These  peaks  occur  in  regions 
where  the  spectral  line  density  of  rotational  lines  of  many  different 
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Fig.  22.  Some  predicted  spectra  at  the  detector  for  a  range  of  atmospheric 
C02  concentrations.  Experimental  C.O^  pressure  is  0.  1  Torr. 
The  standard  atmospheric  mixing  ratio  is  330  ppm. 
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vibrational  bands  is  great  and  significant  self  absorption  is  occurring  (com¬ 
pare  with  Fig.  17).  A  possible  explanation  for  the  over  prediction  of  the 
spectral  intensity  in  these  regions  might  be  the  breakdown  of  the  non¬ 
overlapping  line  approximation  used  in  calculating  the  Voigt  profile  trans¬ 
mission.  If  significant  overlap  were  occurring  the  inclusion  of  non-resonant 
absorption  would  result  in  a  reduction  of  the  predicted  transmitted  intensity 
in  these  regions. 

Experimental  data  taken  with  nominally  0.  2  torr  CO^  present  in  the 
chamber  is  compared  in  Fig.  24  with  a  prediction  using  0.  2  torr  CO^ 
chamber  pressure  and  300  ppm  atmospheric  concentration.  Again,  agree¬ 
ment  is  excellent  except  for  the  intensity  of  the  2350  cm"  *  and  2380  cm  * 
peaks.  The  effect  of  the  purge  under  very  fast  flow  conditions  in  seen  in 
Fig.  25  where  data  from  the  9MAR1  run  is  displayed.  (Note  that  significant 
difference  in  the  spectral  shapes  for  slow  and  fast  purge  data,  see  Fig.  24. ) 
For  comparison,  the  predicted  spectrum  for  the  case  of  a  66  ppm  CO^  con¬ 
centration  in  the  optical  path  is  also  plotted  in  Fig.  25.  The  agreement 
between  the  synthetic  and  experimental  spectra  is  outstanding,  with  only  a 
slight  under  prediction  of  the  intensity  of  the  2310  cm"  *  peak.  Thus,  under 
very  fast  purging  conditions,  a  significant  reduction  in  the  optical  path 
concentration  may  be  achieved.  The  amount  of  COg  present  in  the  chamber 
was  taken  to  be  0.  25  torr  for  the  prediction;  this  adjustment,  in  addition  to 
improving  the  spectral  match,  seems  to  be  in  agreement  with  the  con¬ 

centration  predicted  from  the  CO  relaxation  data.  (See  Section  IV. )  It 
should  be  noted  that  even  under  these  rapid  purging  conditions  the  trans¬ 
mitted  spectral  data  does  not  resemble  the  chamber  window  fluorescence 
(compare  with  Fig.  14  and  15).  Sufficient  absorption  still  occurs  in  the 
atmosphere  to  modify  the  spectrum. 

A  spectrum  from  one  of  the  runs  taken  without  purging  of  the 
optical  path  (5APR2)  is  given  in  Fig.  26.  The  predicted  spectrum  (of 
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Fig.  24  Comparison  of  data  taken  0.  1  msec  before  beam  termina¬ 
tion  with  0.  2  torr  CO2  pressure  (11  MARCH  RUN  1)  and 
prediction  assuming  the  same  chamber  pressure  and 
300  ppm  ambient  CO^  concentration.  Data  is  the  broken 
line,  prediction  is  the  solid  line. 
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Fig.  25  Comparison  of  data,  taken  0.  1  msec  before  beam  termina¬ 
tion  with  0.  1  torr  CO2  pressure  and  a  very  fast  purge,  and 
a  prediction  assuming  0.  25  torr  CO£  chamber  pressure  and 
66  ppm  ambient  CO2  concentration.  Data  is  broken  line, 
prediction  is  solid  line. 
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0.  1  torr  CO^  chamber  pressure  fluorescence)  was  found  to  best  match  the 
data  when  the  atmospheric  CO^  concentration  was  about  1.  3  times  the 
standard  mixing  ratio,  i,  e,  ,  430  +  120  ppm.  The  predicted  spectrum  is 
presented  as  the  solid  curve  in  Fig.  26.  Again,  the  frequencies  of  all 
relevant  spectral  features  are  matched,  although  the  intensity  of  the 
2325  cm-*  shoulder  is  underestimated.  In  general  it  is  concluded  that 
the  detailed  band  signatures  developed  from  the  absorption  modeling  pro¬ 
vide  good  agreement  with  the  data  base. 

It  is  of  course  dangerous  to  rely  solely  on  the  radiation  signature 
comparisons  for  determining  the  ambient  CO.,  concentrations.  Therefore 
a  consistency  check  was  performed  using  the  blackbody  transmission  spectra 
taken  daily  with  the  data.  To  perform*  these  calibration  runs,  the  entire 
interferometer  had  to  be  rotated  so  that  the  blackbody  radiator  could  be 
positioned  at  a  distance  equal  to  that  between  the  chamber  window  and  the 
entrance  aperture  of  the  interferometer  (equalizing  absorption  path  lengths). 
While  rotating  the  interferometer,  the  end  of  the  purged  path  was  open  to 
the  room  air,  and  additionally,  due  to  the  blackbody  aperturing  system, 
there  was  a  5  cm  distance  at  the  blackbody  which  is  unpurged.  As  a  result, 
for  runs  with  a  purging  flow,  the  amount  of  CO^  present  in  the  optical  path 
during  the  calibration  run  is  felt  to  be  at  best  intermediate  between  the 
true  purged  and  atmospheric  values.  The  absorption  feature  in  the 
blackbody  spectrum  is  shown  in  Fig.  27.  There  was  no  attempt  at  optical 
path  purging  for  thus  run. 

In  order  to  generate  a  synthetic  absorption  feature  to  compare  with 
the  calibration,  a  calculation  of  the  absorption  of  radiation  by  each  rota¬ 
tional  line  of  the  three  strongest  CO^  bands  of  Table  I  was  performed  at 

frequency  intervals  smaller  than  the  atmospheric  CO_  linewidth.  Specifi- 

£  _  1 

cally,  at  one  atmosphere,  the  collision  broadened  CO_  HWHM  is  0.  07  cm 

-1  C 

while  the  sampling  interval  was  taken  as  0.  05  cm  .  (It  may  be  seen  from 
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Fig.  27  COz  absorption  feature  in  the  spectrum  of  a  1173  K  blackbody. 
No  purging  of  the  optical  path  was  attempted. 


Eq,  (34)  that  for  this  sampling  interval  the  maximum  absorption  coefficient, 
k^,  calculated  for  a  given  transition  will  necessarily  be  within  10%  of  the 
line  center  value,  kVQ.  )  This  calculation  is  rigorous  in  that  absorption  far 
into  the  wings  of  the  transitions  as  well  as  the  effect  of  overlapping  lines 
are  taken  into  account. 

The  predicted  transmission  of  a  blackbody  source  was  determined 

for  a  range  of  atmospheric  CO-  mixing  ratios.  These  Hi-tran  type  pre- 

-1 

dictions  were  then  convolved  with  a  10  cm  resolution  sine  slit  function. 

The  predicted  atmospheric  CC^  absorption  features  are  plotted  in  Fig.  28 
for  a  range  of  concentrations.  It  is  seen  that  the  fractional  absorption 
changes  rather  slowly  with  CO^  concentration  and  that  the  predicted  curves 
are  very  similar  in  shape  to  the  data  of  Fig.  27,  with  the  absorption  in¬ 
creasing  and  broadening  as  the  amount  of  CC^  is  increased.  Due  to  spectral 
noise,  an  exact  spectral  match  to  one  of  the  predictions  was  not  possible. 
Thus  it  was  felt  that  comparison  of  the  integrated  areas  of  the  absorptions 
would  provide  a  more  exact  ,  ^termination  of  the  atmospheric  CO-,  concen¬ 
tration.  The  observed  fractional  transmission  of  the  blackbody  radiation 
over  the  2310-2380  cm  *  range  is  68%  (assuming  the  detector  exhibits  a 
flat  response  over  this  spectral  range).  The  predicted  fractional  trans¬ 
missions  for  several  atmospheric  CO^  concentrations  are  listed  in  Table  III 
and  it  can  be  seen  that,  the  amount  of  atmospheric  CC^  present  in  the 
unpurged  optical  path  may  be  estimated  to  be  ~  410  ppm  in  reasonable  agree¬ 
ment  with  the  value  deduced  earlier  from  an  analysis  of  the  unpurged 
002(^3)  emission  signature,  see  Fig.  26.  (It  is  seen  however  that  the 
fractional  transmission  is  only  a  weak  function  of  the  atmospheric  CO^ 
mixing  ratio  and  thus  this  is  not  a  very  sensitive  technique  for  determining 
the  mixing  ratio, ) 

Blackbody  calibrations  taken  with  a  purged  optical  path  exhibited  CC>2 
absorption  features  of  the  same  magnitude  as  those  observed  with  unpurged 
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Fig.  28.  Hi-tran  predictions  of  blackbody-  radiation  wit’  atmospheric  CC^ 
absorption  features  for  a  range  of  CO2  press  :res.  In  order  of  C 
increasing  absorption  the  concentrations  are  165  ppm,  330  ppm, 
495  ppm  and  990  ppm. 


TABLE  HI 


PREDICTED  FRACTIONAL  TRANSMISSIONS  OF 
BLACKBODY  RADIATION 


CO 


2 


Concentration 


495  ppm 
430  ppm 
390  ppm 
360  ppm 
330  ppm 
30  6  ppm 
250  ppm 
165  ppm 


Fractional  Transmission 
0.  64 

0.  67 

0.  69 

0.  70 

0.71 

0.  73 

0.75 

0.  82 
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runs.  This  is  probably  the  result  of  the  diffusion  of  C02  ^rom  th®  room  air 
into  the  optical  path  during  interferometer  rotation.  Thus,  the  success  of 
the  purge  in  removing  CO 2  can  only  be  estimated  from  the  emission  data. 
Considering  the  agreement  between  unpurged  concentrations  deduced  from 
the  emission  signatures  and  from  the  blackbody  absorption  data,  this  latter 
method  of  estimating  the  atmospheric  CO 2  is  felt  to  be  reasonably  accurate. 

Finally,  spectral  modeling  of  the  three  high  pressure  runs  (28MAR2, 
18APR2,  and  18APR3)  was  undertaken  in  order  to  observe  the  effect  of 
increased  collisional  broadening  on  the  predicted  spectra.  The  CO^  fluore¬ 
scence  observed  near  beam  termination  for  the  two  18APRIL  runs  is  dis¬ 
played  in  Figs.  29  and  30.  Both  spectra  are  similar  to  those 'observed  at 
lower  pressures.  The  increasing  intensity  at  low  wavenumbers  is  due  to 
the  onset  of  CO  fluorescence  (v  =  1  R  branch)  which  is  the  dominant  spectral 
radiator  under  these  conditions.  Note  that  in  comparison  to  the  low  pres¬ 
sure  data  the  shoulder  on  the  2310  cm"’*  feature  is  more  pronounced,  making 
the  entire  feature  broader,  expecially  for  the  86  torr  case  (Fig.  30).  The 
optical  path  was  not  purged  during  these  runs. 

Test  chamber  transmission  functions  were  tabulated  for  these  cases 

7 

again  using  Whiting's  approximation  for  the  Voigt  lineshapes  with  HWHM 
linewidths,  b^,  of  0.  0087  cm  ^  and  0.  0128  cm  \  for  the  50  and  86  torr 
cases  respectively.  These  transmission  functions  are  plotted  against  SN  L 

XV 

in  Fig.  31  along  with  the  predicted  transmissions  for  the  Doppler  and  Lorentz 
limits  and  the  Voigt  profile  prediction  for  9  torr  taken  from  Fig.  6.  The 
additional  broadening  at  higher  pressures  produces  increased  intensity  in  the 
wings  of  lines  and  thus  increased  relative  transmission.  As  a  result,  the 
626  fundamental  transition  is  the  dominant  contributor  to  the  intensity.  The 
relative  contributions  of  the  626  fundamental  band  and  the  two  strongest  hot 
bands  and  of  the  636  fundamental  band  are  plotted  in  Fig.  32  as  a  function 
of  the  C02  pressure  for  the  85  torr  argon  case.  Note  that  the  fraction  of 
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Fig.  29  Observed  fluorescence,  0.  1  msec  before  beam  termination, 
from  0.  3  torr  CO2,  49.  7  torr  Ar  (18  APRIL  RUN  2)  with  no 
optical  path  purge. 


Fig.  30.  Experimental  data,  0.  1  msec  before  beam  termination,  from  a 

0.  2  torr  CC^,  86  torr  A r  experimental  mixture  (  18  APRIL  RUN  3). 
There  is  no  purging  of  the  optical  path. 
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Fractional  transmission  as  a  function  of  rotational  line- strength  times  column  density 
for  two  high  pressure  Voigt  calculations.  Previous  transmission  functions  are  plotted 
as  broken  lines  for  comparison. 


FRACTION  OF  TOTAL  TRANSMITTED  INTENSITY  AT  CHAMBER  WINDOW 
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C02  PRESSURE  (torr) 

Fig.  32.  Relative  contributions  of  four  of  the  strongest  bands  to  the  total 
transmitted  intensity  at  the  chamber  window  as  a  function  of  CO_ 
pressure  in  85  Torr  of  argon.  ^ 
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transmitted  radiation  due  to  transitions  other  than  the  626  fundamental  band 
is  much  smaller  at  this  increased  pressure  (compare  with  Fig.  9).  This 
effect  is  evident  from  the  predicted  chamber  window  spectra  shown  in 
Figs.  33  and  34  for  test  chamber  CO.,  pressures  of  0.  33  and  0.  10  torr; 

The  hot-band  and  isotopic  radiation  is  apparent  only  as  a  shoulder  on  the 
low  wavenumber  side  of  the  fundamental  transition  P  branch. 

The  total  predicted  CC^v^)  band  intensity  incident  upon  the  detector, 
with  and  -without  the  effects  of  atmospheric  CO^  absorption  are  plotted  in 
Fig.  35;  The  lower  curve  is  calculated  using  an  atmospheric  CO£  concen¬ 
tration  of  430  ppm.  Similar  predictions  for  9  torr  total  pressure  are  sb'iwn 
for  comparison.  It  can  be  seen  that  the  transmitted  intensities.  For  the 
higher  pressure  case  are  substantially  greater  and  vary  less  with  chamber 
CO_  pressure  than  the  9  torr  predictions. 

The  predicted  spectrum  after  convolution  with  a  10  cm  *  resolution 
sine  slit  function  is  plotted  in  Fig.  36,  This  spectrum  was  generated 
assuming  0.  2  torr  chamber  CO^  pressure,  85  torr  Ar  pressure  and  an 
atmospheric  CO^  concentration  of  430  ppm.  Plotted  for  comparison  is  the 
experimental  data  of  Fig.  30.  The  predicted  spectrum  is  similar  in  shape 
to  the  low  pressure  predictions.  The  predicted  shoulder  intensity  is  again 
less  than  the  observed.  However  the  intensities  of  the  2350  and  2380  cm’* 
peaks  are  in  good  agreement  with  the  data,  A  similar  comparison  between 
the  50  torr  data  and  prediction  is  shown  in  Fig.  37.  This  prediction  is  for 
0.  33  torr  of  CO^  in  50  torr  of  argon  with  an  atmospheric  CO^  concentration 
of  430  ppm.  Again,  all  spectral  features  are  matched,  but  the  intensity  of 
the  shoulder  at  2320  cm  *  is  considerably  underestimated.  Considering 
that  the  non- overlapping  line  approximation  becomes  less  appropriate  with 
increasing  pressure  the  comparison  between  data  and  prediction  is 
reasonable. 
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Fig.  34  Predicted  fluorescence  at  chamber  window  for  Voigt  lines 
broadened  by  85  torr  Ar.  Resolution  is  10  cm'l 
PCO2  =  torr.  Maximum  intensity  =  1.  00  x  10*^5 

W /molecule-str-cm" 
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Fig.  36,  Comparison  of  predicted  spectrum  using  0.  2  torr  CO-,  85  torr 
argon  and  430  ppm  atmospheric  CC>2  with  experimental  data 
(broken  line)  of  Fig.  25.  Resolution  is  10  cm"*. 
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Fig.  37.  Comparison  of  predicted  spectrum  using  0.  33  torr  CO2,  50  torr  Ar 
and  430  ppm  atmospheric  CO2  with  experimental  data  (broken  line) 
of  Fig.  24.  Resolution  is  10  cm”^. 
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In  conclusion,  spectral  modeling  has  been  successful  in  matching 
the  experimental  C02  fluorescence  data  over  a  range  of  experimental  C02 
partial  pressures,  experiment  total  pressures,  and  atmospheric  CO 2  con¬ 
centrations  through  consideration  of  attenuation  of  the  emission  by  Voigt- 
broadened  resonant  absorbers  in  the  chamber  and  by  Lorentz  broadened 
absorbers  in  the  detection  system  optical  path.  The  net  result  of  this  analysis 
is  the  quantity  f,  the  fractional  transmission  of  the  CO^v^)  band  fluorescence 
over  the  optical  path.  This  quantity  is  shown  plotted  in  Fig.  38  as  a  func¬ 
tion  of  CO_  test  chamber  pressure  for  P  =9  torr  and  several  atmospheric 
C02  mixing  ratios  and  may  be  used  in  conjunction  with  Eq.  (15)  to  determine 
the  total  CO^Vg)  concentration  within  the  field- of- view. 


E.  Determination  of  the  Excitation  Efficiency 


Total  CO^fv^)  concentrations  observed  within  the  field  of  view  may 
now  be  deduced  from  the  band  intensity  data  through  use  of  Eq.  (15).  This 
has  been  done  for  six  cases  all  corresponding  to  an  argon  pressure  of  9  torr 
and  beam  voltage  of  30  kV  and  the  results  are  displayed  in  Table  IV.  Note 
that  many  of  the  quantities  tabulated  are  double  valued  because  of  uncertainties 
in  the  C02  partial  pressure.  In  each  case  the  first  value  tabulated  is  for  the 
nominal  measured  C02  pressure  which,  in  the  existing  experimental  con¬ 
figuration,  is  difficult  to  measure  accurately  when  CO 2  represents  only  a 
small  fraction  of  the  total  gas  pressure.  The  second  value  listed  is  for  the 
C02  pressure  deduced  from  the  CO  vibrational  kinetic  analysis  presented  in 
the  next  section.  Fortunately,  as  can  be  seen  from  Fig.  38,  the  quantity  f 
varies  only  weakly  with  CO^  pressure  and  thus  this  uncertainty  will  not 
significantly  effect  the  deduced  CO^v^)  concentration. 


The  band  intensities  listed  in  Table  IV  as  well  as  the  deduced  values 

of  N..  are  specific  to  an  irradiation  time  of  2.  4  msec,  just  prior  to  beam 
3 

termination.  The  global  excitation  rate,  ,  may  be  deduced  from  the 
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Fig.  38.  The  specific  transmittivity  as  a  function  of  experimental  CO 
pressure  for  different  amounts  of  atmospheric  CO,:  66  ppm 
(fast  purge),  300  ppm  (slow  purge)  and  430  ppm  (no  purge). 


CO,(v  )  CONCENTRATIONS  AND  EXCITATION  RATES  WITHIN  THE  FIELD  OF  VIEW 

PA  =9  torr  ;  V  =  30  kV 


CO-{v  )  concentrations  by  evaluating  Eq.  (12)  at  the  appropriate  irradiation 
c  3  _3 

time  and  using  the  observed  decay  time  of  Ty^  =  1.  7  x  10  msec.  The 

results  of  this  evaluation  are  also  tabulated  in  Table  IV  and  are  found  to 

1 8 

vary  between  1.  5  -  2.  5  x  10  excitations /sec. 

It  is  of  fundamental  interest  to  compare  this  rate  of  C 0_(v_)  mode 

4  j 

excitation  to  the  total  rate  of  electron  beam  energy  deposition  within  the 
field  of  view.  For  a  thin  target  this  latter  quantity  may  be  defined  as 

PFOV-Ip(i^)Lf  S’  <37> 

-  5  3 

where  p  is  gas  density,  2x10  gm/cm  ,  is  the  diameter  of  the  field 
of  view,  6.4  cm,  and  g*  is  the  fraction  of  beam  energy  deposited  over  the 
distance  which  falls  within  the  circular  field  of  view.  The  quantity  g’ 
has  been  evaluated  to  be  0.  773  for  the  appropriate  beam  shape,  see  Sec¬ 
tion  II.  The  beam  stopping  power,  1/p  d£/dx,  has  been  tabulated  by  Berger 
12  2 

and  Seltzer  as  6.  5  MeV  -  cm  /gm  for  30  keV  electrons  in  Ar.  However, 
it  has  long  been  recognized  that  the  results  of  Be  rger  and  Seltzer  should  be 
considered  as  lower  bounds  because  they  do  not  include  the  effects  of  multi¬ 
ple  scattering.  A  more  accurate  representation  is  given  by  the  Monte  Carlo 
electron  energy  deposition  calculation  performed  using  the  ELTRAN  com¬ 
puter  code.  The  beam  stopping  power  at  the  center  of  the  field  of  view 

,  2 
predicted  from  use  of  this  code  was  14,  5  MeV  -  cm  /gm.  The  beam  stop¬ 
ping  power  does  vary  across  the  field  of  view  and  in  integrated  value  would 
present  a  more  accurate  representation.  This  variation  will  be  included 
in  future  analyses. 

The  rate  of  beam  energy  deposition  within  the  field  of  view  for  the 
six  cases  of  Table  IV  are  tabulated  in  Table  V.  These  listed  values  are 


TABLE  V 


ELECTRON  EXCITATION  EFFICIENCY  FOR  CO-,(v3)  MODE 


PA  =9  torr  ;  v  =  30  kV 
Ar 


Case  ID 

I 

mA 

P  * 

FOV 

watts 

p-3 

watts 

Excitation 
Efficiency  (%) 

9Marl 

0.  59 

.  85 

. 068-. 092 

8.  0-10. 8 

UMarl 

0.  67 

.  96 

.  101-. 107 

10.  5-11.  1 

llMar2 

0.81 

1.  16 

. 104-. 110 

9.  0-9.  5 

l6Mar2 

0.  84 

1.  21 

. 102-. 117 

8.4-9.  7 

l6Mar3 

0.78 

1.  12 

.094-. 107 

8.4-9.  6 

l6Mar4 

0.79 

1.  13 

.092-.  102 

o 

O' 

H 

CO 

Average 

8.  7-10.  0 

Calculated  using  the  beam  stopping  power  of  Ref.  13. 

♦♦First  value  is  for  nominal  C02  pressure,  second  value  is  that 
deduced  from  CO  kinetic  analysis  presented  in  Section  IV. 


specific  for  the  Monte  Carlo  code  value  of  the  beam  stopping  power  at  the 
center  of  the  field-of-view.  Also  listed  in  Table  V  are  the  beam  currents 
for  each  case.  The  total  rate,  within  the  field-of-view,  of  feeding  electron 
energy  to  the  CO^V^)  mode  is  defined  as 


P  =  R  (hv) 
V,  V, 


(38) 


where  hv  is  the  energy  of  a  quanta.  The  total  electron  excitation  effir- 
ciency  for  the  CO^v^)  mode  is  then  just  the  ratio  P^/Pp^y. 

The  resulting  values  of  Pv  and  the  electron  excitation  efficiency 
are  also  shown  in  Table  V  for  the  excitation  rates  listed  in  Table  IV.  The 
deduced  efficiencies  do  not  vary  significantly  from  case  to  case  and  average 
between  8.7-10%  depending  on  the  choice  of  CO_  pressures.  Note  that 
the  fast  purge  case  (9Marl),  which  has  a  band  transmission  correction  f 
about  a  factor  of  two  larger  than  those  of  the  other  cases,  exhibits  an 
efficiency  in  good  agreement  with  the  other  data.  Furthermore  there  does 
not  seem  to  be  any  significant  difference  in  electron  excitation  efficiency 
over  the  factor  of  two  variation  in  CC>2  pressure.  A  large  error  bar  of  +  4% 
must  be  associated  with  these  efficiencies  due  to  uncertainties  in  the  beam 
energy  deposited  and  instrument  misalignment  errors.  The  effects  of  inter¬ 
ferometer  misalignment  on  absolute  intensity  require  further  definition. 


There  does  not  appear  to  be  any  other  measurements  to  which  these 
COz  excitation  efficiencies  can  be  compared.  There  is  however  a  computer 
model  available  which  was  developed  for  determining  the  relative  allocation 
of  electron  energy  deposited  in  gas  mixtures.  This  code  uses  a  set  of  para¬ 
meterized  cross  sections  for  each  gas  to  determine  the  apportionment  of 
absorbed  electron  energy  among  the  various  possible  paths  of  atomic  and 
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molecular  excitations  and  ionizations.  The  calculations  are  performed 

using  the  continuous  slowing  down  approximation  and  include  the  effect  of 

the  total  decay  (to  thermalization)  in  primary  electron  energy  as  well  as 

the  effect  of  all  additional  excitations /ionizations  caused  by  secondary  and 

14 

higher  order  electrons.  This  computer  model  is  well  documented  and 
has  been  previously  used  to  provide  electron  energy  allocation  predictions 
in  pure  Ar^  and  pure  among  other  gases.  Furthermore  the  para- 

17 

meterized  cross  sections  required  as  input  to  the  model  are  now  available 
for  nine  gases:  Ar,  H^,  H2Of  O 2,  N2,  O,  CO,  CO.,  and  He. 

Energy  allocation  predictions  have  been  made  with  this  code  for 
various  mixtures  of  CO^/Ar;  typical  results  are  shown  in  Figs.  39  -  40 
for  the  case  of  1%  C02  *n  Ar,  The  predicted  excitation  efficiencies  for 
various  inelastic  processes  in  Ar  are  shown  in  Fig.  39  as  a  function  of 
primary  electron  energy.  Similar  results  for  C02  are  presented  in  Fig.  40. 
(Notation  is  as  in  Refs.  15  and  16.  )  These  excitation  efficiencies  which  are 
the  fraction  of  primary  electron  energy  partitioned  into  all  relevant  ioniza¬ 
tion,  electronic  state  and  vibrational  excitation  processes,  are  seen  to 
asymptote  for  primary  electron  energies  greater  than  one  keV  and  thus  the 
efficiencies  for  30  keV  primaries  will  be  the  same  as  those  shown  for 
10  keV  primaries. 

The  predicted  excitation  efficiency  for  the  CO^v^)  mode  is  found 
to  be  17.  1%  which  is  nearly  twice  as  large  as  the  measured  values  of  ~  9% 
listed  in  Table  V.  Furthermore  the  predicted  efficiency  is  found  to  vary 
only  slightly  with  C02  mole  fraction,  being  16,  8%  for  a  10%  CO2/90%  Ar 
and  14.  1%  in  pure  CO^.  As  will  be  discussed  in  Section  IV,  the  CO  fluor¬ 
escence  that  is  observed  may  be  due  to  the  dissociative  recombination  reaction 

C02+  +  e  -  CO(V)  +  O  .  (39) 

The  beam  created  species  CO  and  O  remain  in  the  field  of  view  for  several 
pulses,  and  should  be  included  in  the  gas  mixture  composition  for  energy  allo¬ 
cation  predictions.  Inclusion  of  the se  species  decreases  the  predicted  C02 
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(V^)  excitation  efficiencies.  For  example,  at  concentrations  of  0.  03  torr  beam 
created  CO  and  O  would  have  the  effect  of  reducing  the  predicted  efficiency 
to  the  observed  ~9%  value.  At  present  there  is  no  way  to  quantitatively  esti¬ 
mate  the  concentration  of  beam  created  species.  Additionally,  the  CO  exci¬ 
tation  cross-sections  used  in  the  computer  predictions  are  incomplete.  Thus 
the  theoretical  predictions  and  experimental  observations  may  be  reconciled, 
but  accurate  verification  will  require  future  analysis. 

The  CO  pressure  has  only  been  varied  slightly  in  the  present  study, 

^  2 

however  in  previous  work  the  fluorescent  intensity  was  measured  over  a 
factor  of  ten  range  in  CC>2  concentration.  As  an  example  it  was  observed  in 
the  earlier  work  (see  Fig.  60  of  Ref.  2)  that  the  relative  CO^fv^)  band  inten¬ 
sity  per  unit  beam  current  decreased  by  a  factor  of  0.  55  as  the  CO^  pressure 
was  increased  from  0.  1  to  0.  5  torr.  From  Fig.  38  it  may  be  seen  that  the 
predicted  band  transmission  function  f  also  decreases  approximately  by  a  fac¬ 
tor  of  0.  55  over  that  pressure  range.  Thus  the  data,  in  agreement  with  pre¬ 
diction,  are  consistent  with  a  relatively  constant  CO^Vj)  excitation  efficiency 
over  at  least  a  factor  of  five  range  in  C02  pressure.  Additional  runs  over  a 
wide  CC>2  pressure  range  should  be  performed  to  substantiate  this  observation. 


IV.  CO  DATA  ANALYSIS 


The  dominant  infrared  fluorescent  feature  observed  in  electron 
irradiated  mixtures  of  CO^/Ai  is  the  fundamental  vibration/rotation  band 
of  CO.  Fluorescence  from  at  least  the  first  fourteen  vibrational  levels  of 
CO  is  observed  and  the  vibrational  distributions  are  highly  non-Boltzmann 
in  nature.  The  kinetic  interpretation  of  this  data  has  been  considered  in 
great  detail  in  Ref.  2  and  will  only  be  discussed  briefly  here. 

CO  is  not  present  in  the  test  gas  and  thus  must  be  formed  chemically 
through  reactions  involving  electron  beam  created  species.  One  such  pos¬ 
sible  source  would  be  to  the  dissociative  recombination  reaction 

CO*  +  e  -+  CO(v)  +  O  (39) 

Once  the  vibrationally  excited  CO  is  formed  it  is  subject  to  various  de¬ 
activation  reactions,  dominantly  vibrational  exchange  with  CO^ 

CO(v)  +  COz  -*  CO(v-l)  +  co2(v3)  ,  (40) 

and  radiative  decay 

A 

1 

CO(v)  *  CO(v-l)  +  hv  (41) 

A  *> 

CO(v)  CO(v-2)  +  hv  .  (42) 

It  is  also  possible  for  the  vibrationally  excited  CO  to  react  with  electron 
beam  created  species  which  can  build  up  in  concentration  within  the  test 
chamber.  The  most  important  of  these  would  be 
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CO(v)  +  CO<0)  -♦  CO(v-l)  +  CO(l) 


(43) 


where  the  vibrationally  cold  CO  is  the  cumulative  result  of  all  beam  pulses 

occurring  during  the  test  gas  residence  time  within  the  field  of  view.  The 

2 

possible  effects  of  reaction  (43)  have  been  examined  in  some  detail  and 
efforts  have  been  made  to  take  data  under  conditions  where  this  process  will 
be  unimportant. 

A  preliminary  analysis  and  interpretation  of  a  relatively  large  CO 
fluorescence  data  base  will  be  presented  in  the  following  pages.  The  ultimate 
goals  of  this  analysis  are:  (a),  the  determination  of  the  vibrational  level 
dependent  rate  constants  for  reaction  (40),  (b),  the  deduction  of  the  rates  for 
CO(v)  production  within  the  field  of  view,  and  (c),  the  specification  of  the 
electron  beam  efficiency  for  producing  vibrationally  excited  CO. 

The  techniques  used  in  the  data  reduction  and  interpretation  have 
been  discussed  previously2  and  will  only  be  briefly  reviewed  here.  These 
techniques  have  been  used  to  analyze  the  large  data  base  accumulated  during 
the  present  reporting  period.  Only  preliminary  analysis  has  been  performed 
on  the  complete  data  base.  More  detailed  work  is  currently  in  progress, 
and  methods  of  increasing  the  reliability  of  the  deduced  populations  of  high 
lying  vibrational  levels  are  being  investigated.  The  eight  experimental 
runs  which  have  a  suitable  combination  of  CO^  pressure,  flow  rate,  and 
signal  level  to  permit  modeling  of  CO  kinetics  are  presented  in  Table  VI 
along  with  the  relevant  experimental  conditions.  As  can  be  seer  several 
of  the  runs  differ  only  in  the  gas  residence  time  in  the  field  of  view.  No 
variation  in  the  kinetic  relaxation  rate  was  observed  between  these  runs, 
suggesting  that  cold  CO  does  not  play  a  significant  role  in  the  relaxation 
of  the  vibrationally  excited  CO. 

Because  the  vibrational  distributions  are  non- Boltzmann,  determina¬ 
tion  of  the  populations  must  be  made  by  comparison  of  the  intensity  data 


EXPERIMENTAL,  CONDITIONS  FOR  RUNS  INCLUDED  IN  CO  ANALYSIS 


with  a  synthesized  spectrum  which  is  the  sum  of  unit  concentration  radia¬ 
tion  predictions  (basis  functions)  for  each  vibrational  transition,  v*v-l. 

2 

This  comparison  is  made  using  previously  developed  non-linear  regression 
analysis  techniques.  The  synthesized  spectrum  at  any  wavenumber,  I^(v), 
is  the  weighted  sum  of  these  fundamental  transition  basis  functions,  §^(v) 


VV)  =  £  nv  VV}  (44) 

where  the  quantity  n  is  the  total  CO  population  in  level  v  in  the  detector 

field  of  view,  and  v  is  the  spectral  wavenumber.  To  properly  match  the 

spectral  data,  the  basis  functions  must  be  convolved  over  a  slit  function 

representative  of  the  optical  detection  system.  Minimization  of  the  squares 

of  the  difference  between  the  data,  I(v),  and  the  synthetic  spectrum,  i.  e.  , 

2 

(Ip(v)  -  I(v))  ,  permits  determination  of  the  quantities  n^  once  the  basis 
functions  are  specified. 

The  basis  functions  are  comprised  of  the  individual  vibrational/ 
rotational  transitions  of  the  ground  state  CO  molecule.  The  intensities  of 
these  lines  are  given  by  the  relationship 


v->v-  1 
J-*  J+l 


kCVv->v-l  ,  J-»J±1  ^ J -*J±  1 


(2J+1)  v 


3 

v*v- 1 


A 

v-»v- 


1 


(45) 


where  v  is  fixed  for  each  basis  function  and  J  ranges  over  all  thermally 
populated  levels.  In  Eq.  (45)  h  is  Planck's  constant,  c  the  speed  of  light, 

ii  12. 

S„  _  ,  is  the  Honl- London  factor,  A  ,  is  the  Einstein  coefficient,  J 
J-»J±1  ’  v-*v-l 

is  the  rotation  quantum  number,  v  ,  _  _  ,  is  the  line  center  wave- 

n  '  v-*v- 1  ,  J-»J±1 

number  of  the  transition,  V  ,  is  the  wavenumber  of  the  band  origin  and 

V“*V—  1 

N  is  the  population  of  the  radiating  state.  The  quantity  N  may  be 
v,  J  lg  v,  J 

rewritten  as 
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(2Jfl)  exp  (-F  (v,  J)  hc/kT)  , 


(46) 


he  B 


N  _  =  N 


v,  J  v  kT 


where  B  is  the  rotational  constant  for  vibrational  level  v,  defined  as 
v 


B  =  B  -  a  fv  +  1/2) 
v  e  e 


(47) 


and 


F(v,  J) 


J(J  +  1)  B 

v 


D  J2(J  +  l)2 
e 


(48) 


The  line  and  band  origin  wavenumbers  were  calculated  using  the 
anharmonic  oscillator-vibrating  rotor  approximation,  using  the  latest 
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spectroscopic  constants  of  CO  available:  oj  =  2169.812,  uj  x  =  13.289, 

2  e  6  e 
uj  y  =  1.  059  x  10"  ,  B  =  1.93127,  a  =  1.7459  x  10  ,  and  D  =  5.  6268  x 

e  6  -1  e  C  C 
10"  ,  all  in  cm  .  The  Einstein  coefficients  represent  the  least  accurately- 

known  quantity  in  Eq.  (45).  The  values  of  Docken2^  have  been  used  in  the 

present  analysis.  The  intensity  of  a  given  transition  per  vibrationally 

V~*Y—  1 

excited  molecule,  i.  e,  ,  I  /N  may  then  be  found  from  Eqs.  (45)  -  (48). 

J  1  V  | 

Summation  of  these  intensities  over  the  system  slit  function  g(v,  v  )  permits 
construction  of  the  basis  functions: 


5v(v>  =7T  ?  S (v,v 


(49) 


where  the  slit  function  is  as  defined  in  Eq.  (36)  of  Section  IIL  Because 

2 

beam  heating  of  the  gas  is  negligible,  300  K  was  chosen  as  the  tempera¬ 
ture  in  generating  the  basis  functions. 
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Typical  experimental  spectra  from  16MARCH2  (see  Table  VI)  are 
plotted  in  Fig,  41.  The  spectra  are  plotted  at  0.  2  msec  intervals,  and 
cover  the  temporal  range  of  2.  8  -  7.  5  msec,  i.  e. ,  beginning  0.  3  msec  after 
beam  termination.  Note  that  the  spectral  position  of  peak  fluorescent 
intensity  shifts  from  2040  to  1970  cm”*  over  this  time  span;  furthermore, 
the  total  intensity  has  dropped  an  order  of  magnitude.  The  fluorescence 
spectra  occurring  during  the  beam  pulse  are  displayed  in  Fig.  42.  In  this 
figure  time  is  increasing  toward  the  rear  of  the  plot.  During  this  displayed 
time  period,  1.  0  -  2.4  msec,  the  spectral  shape  changes,  with  the  2000  - 
2100  cm”*  region  increasing  in  intensity  more  rapidly  than  the  rest  of  the 
spectrum.  The  rise  in  intensity  above  2220  cm”*  is  due  to  the  onset  of  CO^v^) 
radiation  (see  Section  III).  The  response  of  the  PbSe  detector  used  in  these 
experiments  decreased  rapidly  below  1900  cm”*.  As  a  result  the  response 
corrected  fluorescence  signal  becomes  increasingly  noisy  below  1830  cm  *. 

The  noise  level  in  the  signal  between  1830  and  2200  cm”*  is.  suffi¬ 
ciently  small  to  permit  accurate  fitting  to  the  data  over  that  frequency 
range  for  all  times  displayed  in  the  figures.  The  signal  to  noise  level  in 
this  data  is  significantly  improved  over  that  analyzed  previously.  An 
improvement  in  the  differential  pumping  rate  resulted  in  a  reduction  of  the 
degree  of  beam  spread  prior  to  entrance  into  the  chamber  thus  increasing 
the  electron  beam  current  (and  the  energy  deposited  in  the  field  of  view  ). 
Additionally,  a  solid  state  pulsing  network  was  used  to  improve  beam  pulse 
current  reproducibility.  As  a  result,  the  noise  level  does  not  dominate  the 
signal  until  below  1780  cm”*,  and  evaluation  of  the  populations  of  CO  vibra¬ 
tional  levels  up  to  v  =  14  could  be  attempted. 

Curve  fits  for  all  the  cases  were  generated  over  the  spectral  range 
1780  -  2220  cm"*  and  for  the  time  intervals  shown  in  Table  VI.  The  time 
increment  between  curve  fit  spectra  was  0.  2  msec.  Data  is  available  for 
times  0.  2  -  15.  0  msec  at  0.  1  msec  intervals  for  more  detailed  analysis 
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if  required.  Curve  fits  to  the  data  of  11MARCH2  are  representative  and 
are  shown  in  comparison  to  the  data  in  Figs,  43  -  47  for  different  times 
after  beam  termination.  The  individual  basis  functions  are  plotted  along  the 
bottom  of  each  figure  and  represent  the  relative  contribution  of  each  transi¬ 
tion  to  the  total  intensity.  It  is  seen  that  agreement  is  excellent  not  only 
for  times  near  beam  termination  (Fig,  43)  but  even  at  times  as  long  as 
9.  3  msec  after  beam  termination  (Fig,  47). 

The  spectral  oscillations  in  the  data  are  excellently  matched  in  the 
1900  -  2200  cm”*  region  of  Figs.  43  -  45.  The  noise  level  in  the  signal 
prevents  exact  matching  in  the  1780  -  1900  cm”*  region,  but  the  intensity 
envelope  is  well  matched.  Transitions  with  v=l  1-14  fall  in  this  region,  and 
although  the  calculated  populations  are  reasonably  accurate,  the  estimated 
errors  are  larger  than  for  the  lower,  better  fit  levels.  At  long  times 
(Figs.  46  and  47)  the  intensity  has  decreased  sufficiently  so  that  some  noise 
is  evident  even  in  the  1900  -  2200  cm"*  range,  degrading  the  quality  of  the 
fit  slightly.  The  excellent  agreement  of  the  fits  and  data  in  the  region  of  the 
non- overlapped  v  =  1  R  branch  transition  tends  to  substantiate  the  choice 
of  300°K  as  the  Boltzmann  rotational  temperature. 

Attempts  at  fitting  the  data  down  to  1700  cm”*  yielded  poorer  results 
due  to  non-convergent  oscillations  in  the  predicted  populations.  These  arise 
due  to  the  overlapping  of  the  basis  functions  and  are  an  artifact  of  the  non¬ 
linear  regression  analysis.  The  effect  of  these  large  errors  in  the  deduced 
populations  ofv  =  15  -  17  propagates  back  to  lower  vibrational  states.  The 
present  preliminary  analysis  represents  a  compromise  between  maximum 
use  of  spectral  information  and  minimum  error  propagation,  and  n^  is  the 
maximum  value  used  in  all  subsequent  analysis.  Improved  fitting  methods 
are  currently  being  investigated  to  permit  an  increase  in  the  reliability  of 
the  deduced  populations  of  v  >  14.  The  value  of  n^  is  somewhat  uncertain, 
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Fig.  44.  Curve  fit  to  data  of  11MARCH  RUN  2  at  2.  5  msec  after 

beam  termination  conditions  as  in  Fig.  43.  The  maximum 
intensity  of  the  data  is  1.  54  x  10"  W/str-cm"*, 
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Fig.  45.  Curve  fit  to  data  of  11  MARCH  RUN  2  at  4.  7  msec  after 

beam  termination  conditions  as  in.Fig.  43.  The  maximum 
intensity  of  the  data  is  1.22x10”  W/str-cm”  . 
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Fig.  46.  Curve  fit  to  data  of  11MARCH  RUN  2  at  7.  9  msec  after 

beam  termination  conditions  as  in  Fig.  43.  The  maximum 
intensity  of  the  data  is  1,1  x  10"  W/str-cm"1. 
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Fig.  47.  Curve  fit  to  data  of  11  MARCH  RUN  2  at  9.  3  msec  after 

beam  termination  conditions  as  in  Fig.  43.  The  maximum 
intensity  of  the  data  is  7.  7  x  10”^  W/str-cm“*. 


in  the  present  state  of  the  analysis,  but  small  in  magnitude.  It  is  believed 
that  higher  lying  levels  would  have  even  smaller  populations  and  thus  may 
be  neglected. 

2 

As  observed  previously,  both  the  highest  and  lowest  vibrational 
levels  decay  more  rapidly  than  the  middle  levels.  This  is  graphically 
evident  in  the  intensity  data  as  may  be  seen  by  comparison  of  Figs.  43  and 
47.  The  relative  decays  are  more  clearly  illustrated  in  Fig.  48  where  the 
calculated  populations  of  the  levels  are  plotted  at  various  times  for  the 
11MARCH2  case.  At  the  earliest  time  plotted,  2.  8  msec  (cf  Fig.  43),  the 
population  density  monotonically  decreases  with  increasing  v.  With  increas¬ 
ing  time  the  population  of  level  7  decays  more  slowly  than  the  surrounding 
levels,  and  a  relative  population  inversion  over  levels  4  -  6  is  observed  at 
the  longest  times  shown. 

The  calculated  population  histories  as  a  function  of  time  are  shown 
in  Figs.  49  and  50  for  the  16MARCH3  case.  The  populations  are  plotted  at 
0.  8  msec  intervals  (every  fourth  fit).  The  first  five  levels  decay  linearly 
with  time,  indicating  a  first  order  exponential  decay  rate.  The  broken  line 
represents  the  calculated  population  of  level  1  of  the  data  of  11MARCH1. 

This  data  was  taken  with  a  nominal  experimental  CO ^  pressure  of  0.  2  torr 
and  the  CO  fluorescence  decay  of  this  case  is  observed  to  be  approximately 
twice  as  fast  as  the  16MARCH  RUN  3  data  which  was  taken  with  0.  1  torr 
nominal  CO^  pressure.  Note  that  levels  4-7  are  found  to  decrease  more 
slowly  with  time  than  the  lower  vibrational  levels  and  that  n^  is  approxi¬ 
mately  equal  to  n_  at  long  times.  In  Fig.  50  it  is  seen  that  levels  8-10 

D 

decay  very  slowly.  Furthermore  competition  between  loss  to  lower  levels 
and  feeding  by  the  higher  levels  is  evident;  the  population  of  level  9  actually 
increases  at  times  below  4.  2  msec  as  a  result  of  feeding  from  higher  levels. 
The  highest  levels  shown,  11  -  14,  are  not  fed  by  quanta  from  still  higher 
levels  and  thus  slowly  decay  (although  population  fluctuations  due  to  spectral 
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Fig.  49.  Populations  of  CO  in  levels  v  a  1  -  7  within  the  field  of  view 
as  a  function  of  time  after  beam  onset  for  16  MARCH  RUN 
3  case.  Beam  termination  occurs  at  2,  5  msec.  Broken  line 
is  population  of  level  1  from  11  MARCH  RUN  1  which  has  a 
greater  CO^  density. 
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Fig.  50.  Populations  of  CO  in  levels  v  a  8  -  14  within  the  field  of 

view  as  a  function  of  time  after  beam  onset  for  16  MARCH 
RUN  3  case.  Beam  termination  occurs  at  2.  5  msec. 
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noise  and  fitting  oscillations  tend  to  obscure  the  decay).  The  value  of  n 
plotted  is  artificially  high  at  the  expense  of  level  12.  Additional  data 
scrubbing  and  fitting  will  be  pursued  in  an  effort  to  determine  the  popula¬ 
tion  decay  more  accurately.  It  is  clear  from  Figs.  49  and  50  that  the  very 
slow  decay  rate  of  levels  8-10  provides  for  a  bottleneck  preventing  the 
efficient  relaxation  of  CO  vibrationally  excited  to  v>  8. 

Akinetic  analysis  of  the  relaxation  phenomena  was  performed  using 
the  preliminary  results  for  all  the  runs  of  Table  VL  The  effects  of  CO  self¬ 
relaxation  should  be  small  for  the  predicted  densities  of  beam  created 

2 

CO  present  in  the  field  of  view,  A  more  accurate  estimate  of  these  effects 
requires  a  larger  data  base  covering  a  range  of  beam  created  CO  concentra¬ 
tions.  Uncertainties  in  the  deduced  n^  values  prevent  application  of  the  stan¬ 
dardtime  derivative  rate  equation.  Under  these  circumstances,  the  relaxa¬ 
tion  of  the  vibrationally  excited  CO  by  CO  is  more  accurately  determined  by 

2 

a  time  integral  form  developed  previously: 


n  (t  )  -  n  (t, )  =  a  ,  R  ,  -  a  R  +  A  _  a  _  -  A  a 

v  2  v  1  v+1  v+1  v  v  v+2-*v  v+2  v-»v-2  v 


where  n^t)  is  the  deduced  population  of  level  v  at  time  t  and 


InEq.  (50) 


R  =  k  [CO,]  +  A  , 
v  v  2  v-»v- 1 
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where  k  is  the  rate  constant  for  vibrational  relaxation  by  reaction  (40),  [CO^] 

«•  3 

is  the  CO_  number  density,  cm  ,  and  the  A  .  and  A  _  are  the  funda- 

2  }>  >  v-»v-l  v-*v-2 

mental  and  overtone  transition  Einstein  coefficients  respectively. 

The  direct  determination  of  a^  from  the  vibrational  population  histories 
allows  evaluation  of  Eq.  (50)  for  R^.  Specifically  a  decrease  in  the  total 
population  of  vibrationally  excited  CO  can  result  only  from  fundamental  radia¬ 
tion  from  level  1,  overtone  radiation  from  level  2  and  collisional  relaxation 
from  level  1,  therefore 


yi  -  Vtl))  =  ■  aiRl  “  a2  A2-»0 

v=l 


(53) 


The  remaining  R^'s  can  then  be  found  through  the  relationship 
n-1 

E(n  (t.)  -  n  (t, ))  =  -  a  R  -  A_  a  +  a  R  +  A  a  +A  a 

v  2  v  1  11  2-*0  2  n  n  n+l->n-l  n+1  n+2-+n  n+2 

v=l 

(54) 


Unfortunately,  the  relaxation  data  is  available  only  for  n  <  14,  and  higher 
levels  must  be  neglected  in  the  summations  of  Eqs.  (53)  -  (54).  This  omis¬ 
sion  will  most  seriously  affect  the  upper  levels,  which  were  already  noisy, 
and  levels  8-10  where  the  quantity  (n.^^)  -  n^jt^))  *s  sma^*  Variation  of 
the  number  of  levels  included  in  the  kinetic  analysis  revealed  that  the  solu¬ 
tions  became  more  stable  as  the  number  of  levels  included  in  the  analysis 
was  increased,  as  long  as  the  populations  of  the  higher  levels  were  reasonably 
estimated.  Inclusion  of  populations  which  exhibit  oscillatory  fluctuations 
between  adjacent  levels  or  neglecting  levels  which  have  significant  popula¬ 
tions  such  as  n^,  n^»  or  n  ^  lead  to  physically  unrealizable  values  of  Ry 


(R  <A  ,). 
v  v-»v-l 
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The  values  of  R  determined  for  the  8  runs  of  Table  VI  are  shown  in 
v 

Fig.  51.  Shown  for  comparison  are  the  Einstein  coefficients  The 

R^'s  were  evaluated  over  the  maximum  time  interval  for  which  population 
information  was  available;  this  choice  of  integration  time  minimizes  the 
effects  of  uncertainties  in  the  data.  The  spread  in  the  R  values  is  due  in 
part  to  the  variation  in  experimental  CO^  concentrations  between  cases. 

It  can  be  seen  that  the  relaxation  of  levels  9  -  13  is  dominated  by  the  radia¬ 
tive  emission  term  in  Eq.  (52).  Consequently,  the  relaxation  rates,  k^, 
for  these  levels,  as  determined  from  the  small  differences  between  the 

values  of  R  and  A  . ,  will  have  large  uncertainties  associated  with  them, 
v  v-*v-l 

The  nominal  experimental  CO£  concentrations  present  in  the  irra¬ 
diated  mixture  are  the  values  listed  in  the  second  column  of  Table  VL 

21-2^ 

However,  because  the  relaxation  rate  constant  for  level  1  is  well  known, 
the  actual  CO^  concentration  in  the  field  of  view  may  be  found  by  requiring 
that  the  deduced  kj  equal  the  literature  value.  The  CO£  concentrations 
determined  in  this  manner  are  listed  in  the  last  column  of  Table  VI.  The 
k  values,  calculated  from  R  using  these  latter  values  for  the  CO  number 

V  V  c 

density,  are  plotted  in  Fig.  52.  The  rates  are  smoothly  varying  and  reason¬ 
ably  consistent  between  runs.  The  maximum  relaxation  rates  occur  at  levels 
2  and  3.  The  scatter  in  the  data  increases  with  vibrational  level  beyond 
level  7  reflecting  the  uncertainty  in  the  populations  of  the  higher  levels. 

The  rates  deduced  for  levels  11  -  13  are  essentially  order  of  magnitude 

determinations.  The  broken  line  plotted  in  Fig.  52  represents  the  data  of 
25 

Hancock  and  Smith.  At  large  v  their  rate  constants  fall  at  upper  end  of 
the  scatter  in  the  present  determinations.  The  factor  of  three  difference  in 
the  rate  constants  for  levels  4-6  was  observed  in  our  previous  analysis  and 
is  still  not  understood.  The  present  results  exhibit  good  run-to-run  repro¬ 
ducibility,  however  possible  relaxation  effects  due  to  cold  CO  cannot  be 
completely  precluded. 

These  deduced  relaxation  rates  may  now  be  used  in  the  determination 
of  the  electron  beam  production  rates  of  vibrationally  excited  CO.  The 
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Fig.  52.  Relaxation  rates  for  the  eight  runs  analyzed.  The  experimental 
COg  value  is  chosen  such  that  kj  agrees  -with  the  literature  yalu 
The  broken  line  represents  the  data  of  Hancock  and  Smith.  25 
Runs  are  designated  as  in  Fig.  51. 


fluorescence  data  observed  while  the  beam  is  on,  such  as  that  shown  in 
Fig,  42  was  fit  in  the  manner  described  above.  The  spectra  at  early  times, 
immediately  after  beam  onset,  are  quite  noisy  and  the  quality  of  the  curve 
fit  is  poor.  Typical  fits  are  displayed  in  Figs.  53  -  55  for  the  data  of 
11MARCH1  at  times  of  0.  6,  1.  4  and  2.  4  msec.  The  fit  improves  as  the 
signal  level  increases.  The  calculated  relative  vibrational  populations 
deduced  from  the  fits  to  the  data  of  16MARCH3  are  shown  in  Fig.  56.  The 
populations  of  levels  1-5  increase  linearly  while  the  beam  is  on,  with 
populations  decreasing  with  increasing  vibrational  level.  Levels  6-12 
increase  more  slowly  with  the  exception  of  level  9  which  has  a  low  relative 
population  at  t  =  1.0  msec  and  then  begins  to  increase  rapidly.  It  was  shown 
earlier  that  the  population  of  this  level  continues  to  increase  after  beam 
termination;  thus  the  increasing  in  population  of  this  level  appears  to  be 
due  to  vibrational  cascade  from  higher  levels  rather  than  from  direct  beam 
creation. 

Kinetic  analysis  of  the  excitation  data  is  similar  in  form  to  the  relaxa¬ 
tion,  but  a  term  for  the  absolute  rate  of  creation  of  CO  molecules  in  level 
V,  R-e(v),  must  be  included: 

n  (t)  -  n  (t=0)  =  R  ( v)t  +aR-aR+A_  a  ,  _  -  A  ,  a  (55) 
v  v  '  e  v+1  v+1  v  v  v+2-»v  v+2  v-»v-2  v 

Equation  (55)  may  be  evaluated  directly  inasmuch  as  the  quantities  av  are 
found  from  the  n^  histories  and  R^  was  determined  in  the  relaxation  analysis. 
The  integral  is  taken  over  the  range  0.  0  -  2.  4  msec  (n^ftsO)  =  0),  Because 
the  energy  deposition  within  the  field  of  view  has  only  been  modeled  for  a 
30  kV  electron  beam  this  analysis  was  only  performed  for  the  runs  taken 
at  this  voltage.  The  values  of  R  (v)  deduced  from  these  runs  are  presented 
in  Fig.  57.  The  absolute  creation  rate  is  highest  for  v  =  1  and  decreases 
with  increasing  level. 
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Fig.  53.  Curve  fit  of  1IMARCH  RUN  1  data  at  0.  6  msec  after  beam 
initiation.  Fit  is  performed  over  the  spectral  range  1780  - 
2220  cm"*.  The  basis  functions  ar^  displayed  along  the 
bottom  axis.  Resolution  is  10  cm”  • 
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Fig.  57.  The  absolute  excitation  rates  for  the  30  kV  runs  listed  in 
Table  VI  calculated  using  data  from  the  tim*s  interval 
0.  0  -  2.  4  msec. 


Ill 


If  the  observed  excitation  rates  are  linearly  proportional  to  the  energy 
deposition,  plots  of  R  /I  vs  v  for  all  six  runs  should  be  identical.  (This  is 
because  the  beam  voltages  and  gas  pressures  of  all  cases  considered  are 
the  same. )  Such  a  comparison  is  made  in  Fig.  58  where  the  error  bars  repre¬ 
sent  one  standard  deviation  from  the  average  of  the  data  points.  The  standard 
deviations  are  less  than  10%  for  the  first  seven  levels  and  even  the  four 
highest  levels  have  less  than  30%  error  bars. 

The  sum  of  the  R^fv)  values  for  a  run  represents  the  total  CO  crea¬ 
tion  rate  in  the  field  of  view.  These  values  can  be  compared  with  the  total 
ion  pair  production  rate  in  the  field  of  view  which  is  given  by 


where  the  numerator  is  the  energy  deposition  within  the  field  of  view,  which 
was  calculated  in  Section  HI,  e  is  the  electronic  charge,  and  W  is  the  average 
energy  required  to  produced  an  ion  pair,  26.  1  eV  for  argon.  Values  of 
these  quantities  for  the  six  runs  are  listed  in  Table  VIL  It  is  seen  that 
approximately  four  vibrationally  excited  CO  molecules  are  created  per 
ion  pair.  Additionally,  the  average  excitation  of  a  vibrationally  excited  CO 
molecule  is  found  to  be 


EE  R(v) 

=  7090  ±  260  cm  (57) 

or  approximately  3.  3  quanta  per  vibrationally  excited  molecule.  (E  is 
the  vibrational  energy  of  level  v. )  The  numerator  of  expression  (57)  is 
the  total  beam  power  converted  to  CO  vibrational  energy.  This  quantity, 
which  is  in  the  range  0.  07  -  0. 10  watts  for  the  six  runs,  is  also  tabulated 
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58.  The  effective  excitation  rate?  per  unit  energy  deposition 
for  the  six  30  kV  runs.  Error  bars  indicate  standard 
deviation. 


CO  EXCITATION  PROPERTIES 


in  the  Table.  From  knowledge  of  the  total  energy  deposited,  the  average 
efficiency  for  electron  energy  conversion  to  CO  vibration  excitation  is 
found  to  be  8.  4  +  4.  %.  The  large  error  bars  are  again  due  to  uncertainties 
in  the  energy  deposition  and  interferometer  alignment. 

One  of  the  more  intriguing  aspects  of  this  analysis  is  the  deduction 
that  approximately  2,  5  vibrationally  excited  CO  molecules  are  created  per 
ion  pair.  This  value  is  of  course  inversely  proportional  to  the  energy 
deposited  within  the  field  of  view,  which  requires  more  careful  evaluation. 
Nonetheless  it  is  of  some  interest  to  speculate  on  possible  production 

2 

sources  for  CO.  The  dominant  production  mechanism  assumed  previously 
was  dissociative  recombination  of  CO*,  Eq.  (39),  which  would  provide  for 
at  least  one  CO  molecule  per  ion  pair.  It  has  also  been  found  from  the 
energy  allocation  predictions  of  Section  III,  see  Fig.  40,  that  approximately 
0.  6  electronically  excited  Ar  atoms  are  created  per  ion  pair.  Although  only 
a  fraction  of  these  are  metastable  states,  those  electronic  states  which  do 
have  optically  allowed  transitions  to  the  ground  state  will  be  optically  trap¬ 
ped  at  the  experimental  pressures  and  thus  will  either  radiatively  cascade 
to  metastable  states  or  be  quenched  in  collisions  with  C O^.  These  quenching 
collisions  will  most  likely  be  dissociative  as  has  been  shown  in  the  case  of 
argon  meta stables,  ^  i.  e. , 


Ar(3P)  +  CO  -+  Ar( 1 S)  +  CO  +  O 


Thus  reactions  (39)  and  (58)  can  account  for  at  least  1.  6  CO  molecules  per 
ion  pair. 

Furthermore  both  of  these  reactions  are  highly  exo-ergic  and  have 

26-28  3 

been  observed  to  produce  electronically  excited  states  such  as  CO(a  rr). 

This  latter  metastable  state  will  be  dominantly  quenched  by  collisions  with 

CO_  under  the  present  experimental  conditions,  most  likely  via  the  dissocia- 

«  28 
txve  process 
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(59) 


CO(a3TT)  +  C02  -2  CO  +  O 

CO(a3TT)  formation  has  been  observed  in  CO*  recombination.  Wauchop  and 

28  +  “ 

Broida  found  that  55%  of  the  CO.  created  through  Penning  ionization 

3  3 

collisions  recombined  to  form  CO(a  rr).  If  reaction  (58)  created  CO(a  n) 

•with  unit  efficiency,  the  total  CO  creation  rate  would  be  at  most  2.  75  per 

ion  pair.  Note  that  metastable  O(^S)  could  also  be  produced  by  reactions 

29 

(39)  and  (58).  Preliminary  studies  however,  indicate  that  CO  is  not  pro¬ 
duced  in  quenching  collisions  between  O(^S)  and  CO^. 

Thus  out  of  a  theoretical  maximum  of  2.  75  CO  created  per  ion  pair, 
2.  5  are  experimentally  observed  as  being  created  vibrationally  excited. 

It  must  be  emphasized  that  the  CO  data  analysis  is  not  yet  complete  and 
the  results  presented  above  may  be  revised  when  a  more  detailed  analysis 
is  concluded.  On  the  other  hand,  the  deactivation  rate  constants  and  rela¬ 
tive  excitation  rates  deduced  in  the  present  analysis  are  in  good  agreement 

2 

with  those  determined  previously. 
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V.  SUMMARY  AND  CONCLUSIONS 


An  absolute  calibration  has  been  developed  for  interferometric 
infrared  measurements  on  the  LABCEDE  facility.  The  relative  response 
of  the  system  across  the  field  of  view  has  been  mapped  out  and,  to  first 
order,  is  radially  symmetric  about  the  center  of  the  field  of  view  and  may 
be  described  by  a  Gaussian  distribution.  The  net  effect  of  this  non-uniform 
response,  is  that  effectively  only  one-fourth  of  the  irradiation  resulting  from 
a  uniform  source  filling  the  field  of  view  will  be  detected  by  the  system. 

This  fractional  collection  efficiency  can  of  course  be  different  if  the  fluore¬ 
scence  source  itself  is  non-uniform;  a  range  of  values  has  been  provided  for 
anticipated  electron  beam  shapes.  It  is  clear  that  a  significant  improve¬ 
ment  in  signal  to  noise,  without  a  concomitant  decrease  in  spectral  resolu¬ 
tion,  can  be  achieved  by  introducing  collimating  telescopic  optics  between 
source  and  interferometer.  Such  a  system  modification  would  also  simplify 
data  interpretation. 

The  absolute  calibration  has  been  used  to  analyze  CO^  fluorescence 
data  observed  in  electron  irradiated  mixtures  of  CO^/Ar  in  an  effort  to 
determine  the  fraction  of  deposited  beam  energy  transferred  to  the  CO^v^) 
mode.  The  main  complication  in  this  determination  is  the  evaluation  of  the 
absorptive  losses.  These  have  been  determined  for  the  case  of  Voigt 
broadened  radiative  transitions,  undergoing  absorption  by  Voigt  broadened 
lines  in  the  test  chamber,  under  the  assumption  of  a  non  overlapping  line 
model.  Additional  losses  due  to  absorption  by  ambient  CO^  along  the  optical 
path  have  also  been  evaluated.  Predicted  CO^v^)  spectral  signatures  have 
been  compared  with  measurements  for  conditions  covering  a  range  of  test 
chamber  total  and  CO^  partial  pressures,  and  ambient  CO.,  concentrations. 

In  general  agreement  between  prediction  and  data  has  been  found  to  be  very 
good  with  small  differences  most  likely  due  to  breakdown  in  the  non- overlapping 
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line  assumption.  These  spectral  differences  are  not  anticipated  to  have  a 
significant  effect  on  the  total  predicted  band  intensity. 

Through  this  modeling  it  was  found  that  88  -  94%  of  the  chamber 
CO^CVj)  fluorescence  is  re-absorbed.  By  applying  this  correction  to  the 
measured  band  intensities  it  was  found  that  the  efficiency  of  30  kV  electrons 
for  exciting  the  CO^iv^)  made  in  mixtures  of  9  torr  Argon/0.  1-0.  2  torr  CO^ 
was  9%_+  4%.  Perhaps  the  largest  uncertainties  in  determining  this  exci¬ 
tation  efficiency  are  in  specifying  the  total  beam  energy  deposited  within  the 
field  of  view  and  in  gauging  the  effects  of  interferometer  misalignment  on 
absolute  intensities.  A  careful  evaluation  of  this  quantity  will  be  provided 
in  future  work.  Furthermore,  data  taken  at  other  beam  voltage  and  chamber 
pressure  conditions  remains  to  be  reduced.  It  should  also  be  noted  that  the 
presence  of  beam  created  species  may  provide  for  a  lower  CO_(v.)  excita- 

H  J 

tion  efficiency  than  would  be  obtained  in  pure  mixtures  of  Ar/CO^. 

Preliminary  analysis  of  the  CO  fundamental  vibrational  band  fluore¬ 
scence  data  has  provided  a  set  of  rate  constants  for  the  processes 

CO(v)  +  COz  -»  CO(v-l)  +  COz(v3)  (40) 

The  statistical  uncertainty  in  these  rates  becomes  quite  large  for  vibrational 
levels  greater  than  v  =  7  reflecting  in  part  uncertainties  in  the  populations 
of  the  higher  vibradonal  levels.  The  vibrational  level  dependent  production 
rates  for  CO  have  also  been  deduced.  When  normalized  by  the  rate  of  beam 
energy  deposition,  these  quantities  are  found  to  be  remarkably  consistent 
from  run  to  run.  Preliminary  estimates  are  that  approximately  8  +■  4%  of  the 
deposited  beam  energy  is  transferred  to  CO  vibrational  energy  and  that~2.  5 
excited  CO  molecules  are  created  per  ion  pair  with  average  excitation  level 
of  3.  3  vibrational  quanta.  A  more  careful  analysis  will  include  the  effects 
of  the  beam  created  species  CO  and  O  on  both  the  relaxation  and  produc¬ 
tion  rates.  This  analysis  has  continued  under  a  separate  contract. 
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Lastly,  it  should  be  clear  that  the  development  of  the  absolute  cali¬ 
bration  was  a  complicated  procedure,  involving  both  non-uniform  sources 
and  response  functions,  A  computer  code  to  provide  a  more  exact  solu¬ 
tion  of  the  detector's  response  to  the  non-uniform  electron  beam  source 
should  be  developed.  Given  the  usefulness  of  this  calibration  in  determin¬ 
ing  absolute  creation  rates  and  the  beam  excitation  efficiencies,  it  is 
felt  that  the  process  should  be  repeated  carefully  with  both  experimental 
and  theoretical  crosschecks  made  wherever  possible.  After  this  has  been 
accomplished  the  techniques  discussed  in  this  report  could  be  used  to  de¬ 
termine  other  absolute  electron  beam  excitation  efficiencies  such  as  those 
for  forming  vibrationally  excited  NO  and  in  air. 
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ABSTRACT 


The  infrared  fluorescence  of  electron  irradiated  room  air  and  mixtures 
of  N2/Oz  has  been  observed  over  the  spectral  range  of  1600  -  6700  cm"*. 
Dominant  spectral  features  include  the  NO  fundamental  and  first  overtone 
bands  as  well  as  the  C02(v3)  and  N20(v3)  fundamental  bands.  Comparisons 
between  the  spectra  observed  in  room  air  and  laboratory  N2/02  mixtures 
are  provided. 
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1.  INTRODUCTION 


One  of  the  main  topics  of  discussion  at  the  recent  High  Altitude  Effects 

Simulation  Infrared  Data  Review  meeting*  was  the  origin  of  and  data  base  for 

2.7  ym  band  radiation  occurring  in  the  quiescent  and  aurorally  disturbed  upper 

atmosphere.  The  short  wave  infrared  data  taken  on  several  recent  rocket 

flights  were  discussed  and  some  comparisons  with  theoretical  estimates  were 

provided.  These  predictions  were  developed  under  the  assumption  that  the 

dominant  atmospheric  radiation  source  in  the  2.  7  ym  region  is  the  NO  first 

overtone  vibration/rotation  band.  Although  this  thesis  is  generally  accepted 

as  valid  it  has  not  yet  been  verified  experimentally  inasmuch  as  to  date  upper 

atmospheric  fluorescence  measurements  have  been  performed  utilizing 

radiometers  which  allow  little  spectral  resolution. 

2 

It  was  pointed  out  during  the  meeting  that  there  are  several  sets  of 

measurements  of  2.7  ym  radiation  which  are  not  consistent  with  an  NO 

source.  Specifically,  in  some  instances  measurements  of  the  total  columnar 

2.  7  ym  radiation  as  a  function  of  altitude  exhibit  a  behaviour  typical  of  an 

optically  thick  source,  whereas  atmospheric  NO  is  optically  thin.  Furthermore 

it  was  pointed  out  that  the  ratio  of  2.  7  ym  to  5.'3  ym  band  radiation  observed 

3 

in  recent  project  EXCEDE  measurements  was  much  too  large  to  be  explained 

in  terms  of  the  ratio  of  the  NO  first  overtone  to  fundamental  band  radiation. 

4 

Lastly  some  laboratory  measurements  of  electron  irradiated  air  were  shown 

to  exhibit  spectral  features  due  to  a  species  other  than  NO. 

Of  course  it  can  be  argued  that  the  aforementioned  data  were  either 

anomalous  or  affected  by  contaminants.  Indeed  the  laboratory  spectra  of 

Ref.  4  was  later  found  to  be  contaminated  by  HF,  which  was  formed  in  some 
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complex  manner  from  teflon  spacers  positioned  inside  the  test  chamber. 

In  any  event  the  final  resolution  of  this  conflict  must  await  spectrally 
resolved  measurements  of  the  atmospheric  radiation.  Nonetheless  it  can  be 
of  some  interest  to  define  alternate  atmospheric  sources  of  2.  7  [jm  radiation 
which  would  be  consistent  with  the  data  base. 

2 

Such  an  analysis  has  recently  been  performed  and  it  was  concluded 

that  atmospheric  CC>2  is  the  most  likely  additional  source  of  2.7  (jm  radiation. 

It  is  well  known  that  radiation  from  the  CO^  combination  bands  (101  -+  000)  and 

5 

(021  -»  000)  will  fall  in  the  2.  7  um  region.  These  transitions  are  typically 
taken  to  be  l/25th  as  weak  as  the  corresponding  fundamental  band  transitions 
(101,  021  -»  100,  020).  It  was  pointed  out  however  that  if  these  transitions  were 
only  one  fourth  as  weak  as  the  fundamental  transition  they  could  prove  to  be  a 
significant  source  of  2.  7  |jm  radiation  and  indeed  could  provide  radiation 
signatures  similar  to  the  anomalous  observations  discussed  above. 

The  purpose  of  the  present  study  was  to  examine,  under  laboratory 
conditions,  the  spectrally  resolved  fluorescence  in  the  2.7  |_im  region  occurring 
in  electron  irradiated  air,  and  to  identify  the  molecular  sources  of  this 
radiation.  This  was  done  utilizing  the  room  temperature  LABCEDE  facility 
which  requires  total  operating  pressures  greater  than  one  torr.  It  must  be 
emphasized  that  this  device  does  not  simulate  upper  atmospheric  conditions. 
Specifically  a  high  degree  of  collisional  quenching  will  occur  in  the  test 
chamber  and  -.his  quenching  can  obscure  radiative  phenomena  which  might 
be  dominant  at  the  lower  pressures  of  the  upper  atmosphere. 
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II.  EXPERIMENTAL 


The  LABCEDE  facility  is  a  room  temperature  apparatus  designed  to 
study  the  fluorescence  of  electron  beam  irradiated  gas  mixtures.  In  this 
device  a  constant  pressure,  continuously  flowing  gas  is  irradiated  by  a 
collimated,  mono- energetic  electron  beam.  The  infrared  fluorescence  from 
the  irradiated  gas  target  is  monitored  in  a  plane  perpendicular  to  the  excita¬ 
tion  beam.  This  fluorescence  may  be  both  spectrally  and  temporally  resolved 
through  use  of  a  computer  interfaced  Michelson  interferometer.  Because  of 
signal  to  noise  limitations  the  device  cannot  be  operated  at  gas  pressures  below 
approximately  one  torr  and  thus  cannot  be  used  to  provide  an  exact  simulation 
of  processes  occurring  in  the  upper  atmosphere,  e.  g.  H  ?>  80  km.  In  the 
present  studies  the  observation  region  flow  and  pressure  conditions,  as  well 
as  beam  current  and  voltage,  were  chosen  to  maximize  the  fluorescence 
intensity  in  the  short  wavelength  infrared  (SWIR)  near  2.  7  Hm. 

6-9 

The  LABCEDE  apparatus  has  been  described  in  detail  elsewhere,  and 

only  an  overview  of  the  facility  will  be  presented  here.  The  electron  beam  is 

operated  in  a  pulsed  mode  and  in  the  present  study  was  characterized  by  a  square 

wave  pulse  of  2.  5  msec  duration  and  period  (25  msec  10%  duty  cycle).  The 

electrons  are  accelerated  through  an  adjustable  32  -  36  kV  potential  drop. 

focussed,  and  sequentially  passed  through  three  pin  hole  nozzles  which  form  a 

differential  pumping  network.  The  differential  pumping  permits  observation 

chamber  operation  at  pressures  of  up  to  150  torr  while  maintaining  pressures 
_  3 

below  10  torr  in  the  electron  gun.  The  electron  beam  current  reaching  the 

observation  chamber  is  largely  determined  by  the  amount  of  beam  spreading 

prior  to  the  last  pinhole  nozzle.  The  flow  and  pressure  in  the  observation 
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region  directly  affect  the  nozzle  chamber  pressure  due  to  pump  loading,  and 
thus  can  affect  the  amount  of  beam  scattering  loss  for  a  given  run.  Typical 
electron  beam  currents  were  1-2  milliamps  for  32  -  36  kV  operation.  The 
beam  is  maintained  at  approximately  constant  operating  conditions  during  the 
course  of  a  run;  short  term  stability  fluctuations  and  long  term  drift  were  both 
less  than  5%  of  the  total  current  level. 

The  gases  are  introduced  into  the  chamber  above  and  to  the  side  of 
the  beam  entrance.  Mixing  of  the  test  gases  occurs  just  outside  of  the  observa¬ 
tion  chamber.  The  gas  flow  rate  is  measured  with  a  calibrated  flowmeter 
and  accurately  controlled  by  knife-edge  needle  valves.  The  chamber  pressure 
is  monitored  by  MKS  Baratron  capacitance  manometer  with  an  atmospheric 
pressure  head,  referenced  to  the  electron  gun  pressure. 

The  fluorescence  is  monitored  by  a  lead- salt  detector  mounted  at  the 
exit  plane  of  the  scanning  Michelson  interferometer.  The  moving  interferometer 
mirror,  mounted  on  an  oil  driven  piston,'  is  continuously  scanned.  Relative 
mirror  movement  is  monitored  by  the  interference  fringes  resulting  from  a 
helium- neon  laser  beam  positioned  parallel  to  the  fluorescent  signal.  Con¬ 
structive  interference  fringes  occur  every  time  the  mirror  changes  position 
by  0.316  microns.  The  mirror  scanning  rate  is  chosen  so  that  on  the  average 
40  beam  pulses  occur  between  laser  fringes.  The  detector  signal  is  bandpass 
amplified  (PARC  Model  113),  interfaced  with  a  Digital  PDP-15  computer 
using  an  A/D  converter,  processed,  and  stored  on  magnetic  tape. 

Observation  of  the  temporal  dependence  of  the  fluorescence  spectrum 
is  made  possible  by  a  sample  and  hold  network  referenced  to  the  electron 

Q 

beam  pulse  onset.0  The  network  output  is  directly  input  into  the  A/D  converter. 


at  sampled  intervals  less  than  the  detection  system  response  time.  The  A/D 
converter  time  resolution  limitation  is  around  40  ^seconds,  and  for  these 
experiments,  the  resolution  limit  is  determined  by  the  detector  and  amplifi¬ 
cation  network.  Because  of  the  detector  limitations  it  was  not  possible  to  measure 
the  emission  of  the  NO  fundamental  and  first  overtone  vibration/rotation  bands 
simultaneously,  A  PbSe  detector,  cooled  to  77°K,  was  used,  to  monitor 
fluorescence  in  the  mid  wavelength  infrared  region  (MWIR)  near  5  |Jtm.  This 
detector  had  a  75  jjsec  response  time  which  was  degraded  to  200  nsec  by  the 

bandpass  amplifier  in  order  to  decrease  the  noise  level.  A  PbS  detector  (at 
o 

195  K)  was  used  for  the  SWIR  measurements.  This  detector  exhibits  higher 
responsivity  but  is  limited  to  a  2  millisecond  response  time. 

Averaging  of  the  appropriate  detector  intensities  between  laser  fringes 
is  performed  to  increase  signal  to  noise.  These  averaged  signal  are  then 
stored  on  magnetic  tape  as  a  single  point  in  each  of  the  time  resolved  inter- 
ferograms.  For  all  the  experiments  sufficient  data  points  (at  different  mirror 
positions)  were  gathered  to  provide  for  10  cm  ^  spectral  resolution.  Degrada¬ 
tion  of  the  resolution  to  artificially  decrease  the  apparent  noise  level  could  then 
be  performed  mathematically  by  the  computer  at  a  later  time  if  desired. 

The  gases  used  in  these  experiments  were  ultrahigh  purity  nitrogen 
and  oxygen,'  and  room  air.  When  use  was  made  of  room  air,  it  was  introduced 
into  the  chamber  through  the  flowmeter  and  needle  valve  in  order  to  match  the 
cylinder  gas  flow  conditions  as  closely  as  possible.  The  room  air  could  be 
filtered  to  remove  dust  and  water  vapor  by  means  of  a  glass  wool/silica  gel 
desiccant  trap.  This  trapping  system  was  found  to  be  effective  in  removing 
water  vapor  from  the  laboratory  air. 
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The  infrared  fluorescence  intensity  was  found  to  be  optimized  at  relatively 
slow  flow  conditions,  2.  5  standard  liters  per  minute  total  flow,  with  pumping  of 
the  test  chamber  occurring  only  through  the  electron  beam  differential  pumping 
system.  The  optimized  chamber  operating  pressure  was  100  torr  and  the 
residence  time  of  the  test  gas  was  30  seconds.  Under  these  flow  conditions  beam 
created  species  are  likely  to  build  up  in  the  observation  region.  Specifically 
nitrogen  and  oxygen  atoms  are  produced  in  the  test  gas  either  by  direct  electron 
impact  dissociation  or  by  subsequent  electron/ion/neutral  reactions.  These 
atoms  can  then  produce  other  chemical  species  such  as  NO,  O^,  NC>2,  N^O. 

The  present  experiments  are  designed  to  study  NO  fluorescence  resulting  from 
the  reactions 


N(4S) 

+  Oz  -»  NO(v)  +  O 

(1) 

N(2D) 

+  Oz  •*  NO(v)  +  O 

(2) 

If  the  gas  residence  time  is  long  compared  to  the  beam  pulsing  time  the  amount 
of  NO  and  other  beam  created  species  will  increase  pulse  by  pulse  until  reaching 
a  steady  state  level.  The  details  of  the  relevant  chemical  kinetics  for  N2/02 
mixtures  may  be  found  in  Ref,  6,  It  is  not  anticipated  that  the  presence  of  these 
beam  created  species  will  affect  the  conclusions  of  this  study  however  their 
presence  must  be  considered  when  interpreting  the  experimental  observations. 

Attempts  to  decrease  the  gas  residence  time  result  in  a  loss  of  fluore¬ 
scent  intensity.  Increasing  the  flow  and  pumping  rates  at  constant  chamber  pressure 
lead  to  beam  current  degradation;  pressure  reduction  at  constant  flow  causes  a 
reduction  in  the  total  energy  deposition  of  the  beam  within  the  field  of  view. 
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As  a  result,  only  slow  flow  measurements  were  made  in  the  MWIR  region.  On 
the  other  hand,  because  of  the  greater  sensitivity  of  the  PbS  detector,  SWIR 
fluorescence  measurements  could  be  made  both  at  the  slow  flow  conditions  and 
also  at  flow  rates  of  5  standard-liters  per  minute.  For  these  latter  flow  rates, 
and  a  chamber  pressure  of  12  torr,  the  beam  current  degradation  may  be  kept 
at  an  acceptable  level,  while  the  gas  chamber  residence  time  is  decreased  to 
2  seconds.  Consequently,  in  these  measurements  the  build-up  of  beam 
created  species  is  greatly  reduced.  Experiments  performed  in  the  past^  with 
an  integrating  sphere  in  the  chamber  have  resulted  in  increased  signal  level 
at  the  expense  of  the  residence  time  and  field- of- view  definition.  No  attempt 
was  made  to  implement  the  integrating  sphere  in  the  present  study. 

Experimental  measurements  were  first  made  in  the  MWIR  region  in 
order  to  observe  the  temporal  behavior  of  the  NO  fundamental  band  radiation 
(Av  =  1)  with  0.2  millisecond  resolution.  Additionally,1  the  presence  and 
relative  strengths  of  radiation  from  COgCVj)  and  NzO(v3)  transitions  at  4.3  |_im 
and  4.  5  |_Lm  could  be  examined.  A  series  of  slow  flow  runs  using  first  a 
mixture  of  80%  N2/20%  O-^then  trapped  (dried)  room  air,  and  finally  unfiltered 
room  air,  were  made  keeping  chamber  pressure,  beam  current  and  voltage 
constant  at  100  torr,  1,6mA,  and  36  kV  respectively.  When  unfiltered  air 
was  used,  no  significant  detector  signal  could  be  obtained.  It  is  expected 
that  the  level  of  fluorescence  in  this  case  is  diminished  because  of  the  efficient 
quenching  of  vibrationally  excited  states  by  water  vapor. 

A  similar  series  of  measurements  was  made  in  the  SWIR,  NO  overtone, 
spectral  region  with  the  PbS  detector.  The  experimental  conditions  of  the  MWIR 
measurements  were  duplicated  as  much  as  possible.  Beam  currents  achieved 
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were  slightly  different,  1.  2  mA  -1.8  mA,  but  all  other  parameters  were 
unchanged.  Because  of  the  increased  sensitivity  of  this  detector,  signal 
was  observed  even  for  the  case  of  unfiltered  air.  Fast  flow  measurements 
were  made  with  both  cylinder  gases  and  filtered  room  air  at  36  kV,  but  at  a 
reduced  beam  current  of  0.  9  -  1.4  mA.  Although  the  magnitude  of  the  beam 
current  was  decreased  in  these  runs,  the  current  stability  was  maintained 
at  the  5%  fluctuation  level.  The  total  fluorescence  intensity  for  these  cases 
was  reduced  by  more  than  a  factor  of  two  below  that  of  the  slow  flow  cases. 

An  additional  set  of  measurements  at  reduced  oxygen  partial  pressures 

were  made  both  in  the  MWIR  and  SWIR  regions.  At  the  same  total  pressure, 

it  was  found  that  mixtures  with  fractional  oxygen  concentrations  of  1,  2  - 
-4 

3.  0  x  10  of  the  nitrogen  pressure  provided  the  maximum  NO  fluorescence 
intensity.  The  total  fluorescence  intensity  of  these  runs  was  a  factor  of  four 
larger  than  those  observed  in  the  air  mixtures. 

The  final  collected  data  base  included  eight  runs  taken  with  the  PbSe 
detector  and  fifteen  runs  taken  with  the  PbS  detector.  In  several  instances 
measurements  were  repeated  over  a  period  of  several  weeks  in  order  to  ensure 
data  reproducibility.  Furthermore  detector  spectral  response  calibrations  were 
performed  daily  using  a  black  body  radiation  source.  The  collected  experi¬ 
mental  data  (stored  on  magnetic  tape)  was  processed  by  first  sorting  the  data 
into  time-resolved  interferograms,  Fourier  transforming  them  into  spectra, 
then  applying  a  correction  for  detector  spectral  response  and  field-of-view. 

The  final  form  of  the  reduced  data  is  a  series  of  time- re  solved,  calibrated 
spectra  for  each  run.  Spectral  resolution,  up  to  the  maximum  of  10  cm’\  is 
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d.tarmined  by  the  number  of  interferogram  points  included  in  the  transformation 
These  spectra  are  then  used  for  species  identification  and  kinetic  anaiysis. 


in.  DATA  ANALYSIS  AND  INTERPRETATION 


The  signal  to  noise  ratio  of  the  data  taken  in  the  MWIR  was  suf¬ 
ficiently  low  so  that  the  observed  spectra  exhibited  a  "spikey"  structure. 

The  apparent  noise  level  could  be  diminished  by  artificially  degrading 
the  spectral  resolution  to  20  cm  however,  this  had  to  be  done  with 
care  inasmuch  as  real  spectral  features  can  be  inadvertantly  washed  out 
by  this  procedure. 

A  typical  MWIR  spectra,  taken  in  a  mixture  of  80%  N^  /20%  at 

a  total  pressure  of  100  torr,  is  shown  in  Fig.  1.  This  data  was  taken 

with  the  PbSe  detector  under  slow  flow  conditions  and  corresponds  to  a 

time  just  prior  to  beam  termination;  the  spectral  resolution  of  the  data 

has  been  degraded  to  20  cm  *.  Two  distinct  spectral  features  are 

apparent.  The  first  of  these,  spanning  the  spectral  range  of  1760-1925  cm  *, 

is  the  result  of  NO  fundamental  vibrational/ rotational  band  fluorescence, 

arising  from  NO  molecules  with  up  to  six  quanta  of  vibrational  excitation. 

Radiation  from  the  ^O  (v^)  band  is  also  observed  at  2223  cm  *.  Note 

that  neither  of  these  species  is  present  in  the  original  test  gas  and  thus 

both  must  be  created  through  chemical  reactions  such  as  (  1  and  (2)»  From 

5,  10 

a  knowledge  of  the  relevant  Einstein  coefficients  ,  it  can  be  determined 
that  the  populations  of  the  vibrationally  excited  NO  and  N^O  responsible 
for  the  fluorescence  are  in  the  ratio  of  seven  to  one.  No  other  significant 
spectral  radiators  are  observed  in  the  bandpass  of  this  detector.  (The 
feature  at  1630  cm  *  is  an  artifact  of  the  spectral  calibration.) 
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The  same  spectrum  is  shown  in  Fig.  2  at  10  cm  resolution.  The 
large  noise  level  is  due  to  the  combination  of  poor  detector  response  and 
low  fluorescence  levels.  It  can  be  seen,  however,  that  the  basic  spectral 
features  of  Figs.  1  and  2  are  the  same.  The  spectrum  of  Fig.  3  corres¬ 
ponds  to  the  same  case  at  a  time  of  0.2  msec  after  beam  ter,  ination.  The 
time  difference  between  Figs.  2  and  3  is  equal  to  the  response  time  of  the 
detection  network.  The  intensity  scale  is  the  same  for  both  figures  and  it 
can  be  seen  that  with  increased  time  the  fluorescence  due  to  NO  is  greatly 
diminished,  while  the  N^O  signal  is  essentially  unchanged.  The  rate  con¬ 
stants  for  vibrational  quenching  of  NO  (v=l)  by  N^  and  O^  have  been  mea¬ 
sured**  as  1.7  x  10  *^  cm*/sec  and  2,4x  10  *^  cra*/sec  respectively 
leading  to  a  predicted  NO  fluorescence  decay  time  of  0.07  msec  for  the 

present  case.  The  observed  NO  decay  is  limited  by  the  time  resolution  of 

12 

the  experiment.  N^O  (v^)  quenching  by  N^  has  been  measured  to  have  a 

-15  3 

rate  constant  of  4  x  10  cm  /sec,  and  thus  the  N^O  (v^)  fluorescence 

-4 

would  be  expected  to  decay  on  the  same  time  scale,  as  NO,  10  seconds. 
The  N^O  fluorescence  in  this  experiment  is  observed  to  decay  with  a  time 
constant  of  20  msec,  thus  implying  the  presence  of  a  source  of  N  O  vibra- 

M 

tional  excitation  in  the  test  gas  after- beam  termination.  Vib rationally 

excited  molecular  nitrogen  (i.  e, ,  excited  by  beam  electrons)  is  a  known 

6 

energy  reservoir  under  these  conditions  and  near-resonance  vibrational 
exchange  between  ^Ofv^)  and  N^fv)  is  postulated  as  an  explanation  for  the 
observed  slow  N^O  decay.  Specifically,  the  vibrational  exchange  reaction 
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(3) 


N2(v=1)  +  N2<D  «-»  N20(v3)  +  N2  AE  =  137  cm 
is  assumed  to  be  in  local  steady  state  so  that  to  first  order 

N2<Y=1>  N2°<y3>  -AE/kT 

n2  ‘  n2o  6 

and  thus,  since  N20  is  a  trace  species,  the  decay  of  N20(v3)  mirrors  the  decay 
of  vibrationally  excited  N2»  The  large  energy  defect  between  N2  and  NO, 

484  cm-*,  precludes  a  similar  effect  in  the  case  of  NO. 

In  Fig.  4  the  spectrum  of  irradiated  filtered  room  air  under  slow-flow 
conditions  is  displayed  at  beam  termination.  Both  NO  and  ^Ofv^)  fluorescence 
are  observed  along  with  radiation  from  C02(v3)  at  2350  cm-*  (4.  3  |jm)  as 
expected.  The  N20(v3)  fluorescence  intensity  is  increased  by  20%  over  that 
of  the  runs  of  Figs.  2  and  3,  while  the  NO  fluorescence  has  decreased 

by  a  third.  Carbon  dioxide  is  the  dominant  spectral  radiator.  From  considera¬ 
tion  of  the  relative  bandstrengths  of  CO  ,  N_0,  and  NO,  it  may  be  estimated 

Ct  « 

♦ 

that  the  N20(v3)  and  NO  concentrations  at  b  Jam  termination  relative  to  that  of 

^^2^v3^  are  an<*  5,  6  to  1  respectively.  A  spectrum  from  the  same  case 

corresponding  to  the  fluorescence  0.2  msec  after  beam  termination  is  shown 

in  Fig.  5.  Again  the  NO  fluorescence  has  disappeared  while  both  the  N20  and 

C02  fluorescence  remain  unchanged.  The  fluorescence  detected  3.4  msec 

after  beam  termination  is  displayed  in  Fig.  6.  Vibrationally  excited  N20  and 

C02  are  still  both  present  in  the  sam>  relative  concentrations.  The  measured 

decay  time  of  their  fluorescence  is  approximately  14  msec.  Like  N20,  C02 

exhibits  a  rapid  near  resonant  vibrational  exchange  reaction  with  Ng,  i.  e. , 
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intensity  is  9,  1  x  10" 


N2(v=1)  +  C02  C02(v3)  +  N2  AE  =  19  cm-1  (5) 

and  thus  the  CO^  fluorescence  will  also  mirror  the  decay  of  vibrationally 
excited  N2> 

A  spectrum  taken  under  slow  flow  conditions  in  filtered  room  air 
using  the  PbS  detector  is  shown  in  Fig.  7.  This  spectrum,  corresponding 
to  the  time  of  beam  termination,  has  a  resolution  of  20  cm  1  and  is  uncor¬ 
rected  for  the  spectral  response  of  the  detector;  the  CO,  v,  radiation  at 

M  J 

2350  cm  is  by  far  the  strongest  feature  observed  although  this  detector  is 

relatively  insensitive  to  4.3  jam  radiation.  N20(v3>  fluorescence,  not  a 

prominent  feature  in  the  spectrum  shown,  has  also  been  observed  with  the 

PbS  detector.  The  calibrated  intensity  of  the  CO^Vj)  radiation  shown  in  Fig.  7 

agrees  to  within  20%  with  a  value  obtained  using  the  PbSe  detector  under  similar 

experimental  conditions.  This  excellent  agreement  *'eflects  both  the  accuracies 

of  the  detector  calibrations  and  the  good  run  to  run  reproducibility  of  the  data 

base.  The  other  spectral  features  evident  in  Fig.  7  are  the  radiation  of  interest 

in  the  spectral  range  of  3500  -  3800  cm  *  and  the  sharp  features  around  5500  cm' 

(1.8  Hm)  which  have  been  identified  as  atomic  oxygen  and  nitrogen  lines  in 

13 

previous  studies  performed  at  the  Air  Force  Geophysics  Laboratory. 

Figure  8  displays  a  calibrated  version  on  the  spectrum  of  Fig.  7  over 
the  spectral  range  of  3100  -  3800  cm"1  (3.23  -  2.63  |am).  Also  shown  for 
comparison  is  a  computer  generated  NO  first  overtone  spectrum.  The  generated 
spectrum  is  calculated  using  the  latest  spectroscopic  constants  of  NO  and 
includes  contributions  from  both  total  angular  momentum  states.  The  shapes 
of  all  vibrational/rotation  transitions  in  the  synthetic  spectrum  are  broadened 
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by  a  sine  slit  function  appropriate  to  the  interferometer  used  in  these  experi¬ 
ments.  The  synthesized  spectrum  was  chosen  to  have  a  Boltzmann  vibrational 
distribution  corresponding  to  a  vibrational  temperature  of  5000°K  and  a  rota¬ 
tional  temperature  of  300°K;  no  attempt  has  been  made  to  fit  the  two  spectra. 

It  is  seen  that  the  structure  of  the  experimental  spectrum  taken  in  room  air 
lines  up  very  well  with  the  predicted  NO  spectrum,  leading  to  the  conclusion 
that  NO  is  the  dominent  radiator  in  the  2.  6  -  3.2  (_im  region  under  these  con¬ 
ditions.  Again  fluorescence  arising  from  vibrationally  excited  NO  up  to  level 
v=6  is  observed. 

When  the  spectral  resolution  is  increased  to  10  cm  *  reproducible 
features  not  attributable  to  NO  appear,  as  is  shown  in  Fig.  9.  Although  the 
noise  level  of  the  10  cm"*  calibrated  spectrum  is  greater  than  the  20  cm"* 
spectrum  of  Fig.  8,  the  features  indicated  by  arrows  are  reproducible  both  as 
a  function  of  time  and  from  run  to  run.  The  expected  spectral  shapes  of  the 
^0(02°  1  -*  00°0)  and  (10°1  -*  00°0)  combination  bands  as  well  as  the  CC>2(02ol  -» 
00  0)  and  (10°1  -*  00°0)  combination  bands  are  also  shown  in  Fig.  9.  It  is  clear 
that  CO£  radiation  does  not  occur  in  any  significant  amount  in  this  spectral 
region.  The  N^O  combination  band  at  3480  cm  *  does  seem  to  match  the 

frequency,  if  not  the  shape,:  of  some  of  the  unidentified  spectral  features. 

5 

(This  band  is  predicted  to  be  the  strongest  in  this  spectral  region  for  N^O 
in  vibrational  equilibrium  at  room  temperature.)  Nonetheless,  even  if  the 
3480  cm  *  features  are  attributed  to  N.O,  other  strong  features  remain 
unexplained. 

The  SWIR  spectrum  of  a  slow  flow  run  taken  with  a  80%  N.,/  20% 
mixture  is  shown  in  Fig.  10.  The  data,  shown  with  10  cm  *  resolution,  is 
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taken  at  the  time  of  beam  termination.  Again,  as  for  the  filtered  room  air 
cases,  NO  is  the  dominant  radiator  in  the  spectral  region  shown.  The  lack 
of  significant  CO^  radiation  in  this  spectral  region  is  demonstrated  in  Fig.  11 
which  provides  a  comparison  between  the  spectra  recorded  in  filtered  room 
air  and  an  80%  ^,^20%  mixture  where  no  CO^  should  be  present.  It  can 
be  seen  that  there  are  no  observable  differences  between  the  two  scans,  with 
the  exception  that  the  laboratory  air  case  appears  to  have  a  somewhat  greater 
spectral  resolution.  Interestingly  enough  even  the  features  attributed  to  noise 
appear  to  line  up;  the  unexplained  spectral  features  of  the  air  runs  are  also 
present  in  the  N^/O^  experiments. 

Measurements  taken  under  faster  ilow  conditions  and  for  much  shorter 
residence  times  yield  similar  spectra.  In  these  cases  only  runs  using  the 
PbS  detector  were  possible  due  to  decreased  signal  intensity,  but  the  unidentified 
features  remained.  If  these  features  were  attributable  to  a  build  up  of  beam 
created  species,  such  as  N^O,  their  relative  intensity  should  have  been  greatly 
diminished  because  of  the  decreased  residence  time.  In  fact,  ^Ofv^)  radia¬ 
tion  at  4.  5  |im  was  no  longer  observed  for  these  runs  (with  the  PbS  detector). 
Radiation  arising  from  NO  Av=2  transitions  is  still  the  dominant  radiation 
source  at  2.  7  ^m  under  these  conditions,  and  the  NO  vibrational  distribution 
appears  to  be  very  similar  to  that  of  the  higher  pressure  runs. 

The  experiments  performed  using  trace  concentrations  of  oxygen 
mixed  with  nitrogen  provided  slightly  more  signal  in  the  MWIR  region.  No 
evidence  of  N^O  fluorescence  at  4.  5  |im  is  observed  in  these  spectra.  The 
observed  vibrational  distribution  of  NO  is  not  greatly  different  from  that  of 
the  80%  N2/20%  O ^  mixture  data  of  Fig.  2.  In  these  runs  the  nitrogen  pressure 


-157- 


Comparison  of  N_/©_  and  room  air  spectra  taken  under  similar  conditions.  Resolution  is 
10  cm"^.  Lower  trace:  room  air  fluorescence  (Fig.  9).  Upper  trace:  N_/Q  fluorescence 


was  varied  from  25  -  115  Torr  and  the  oxygen  pressure,  set  to  achieve  maximum 
intensity,  was  varied  from  0.  2  -  0.75  Torr.  No  variation  in  the  observed  NO 
excited  vibrational  distribution  was  noted  over  this  range  of  conditions.  The 
measured  NO  fluorescence  decay  rate  was  in  good  agreement  with  previous 
observations.**  A  more  systematic  study  of  this  type  might  permit  the  measure¬ 
ment  of  individual  vibrational  level  decay  rates  of  NO  by  N^  and  O^. 

When  the  PbS,  SWIR  sensitive  detector  was  used,  the  observed  NO 
spectral  intensity  increased  by  a  factor  of  two  over  the  air  cases.  The 
unidentified  features  noted  previously  all  increased  by  an  even  greater  amount, 
about  a  factor  of  four,  and  are  the  most  prominent  spectral  features,  as  shown 
in  Fig.  12.  The  displayed  spectrum  is  from  a  100  torr  N£,  0.  03  torr  O^ 
slow  flow  run,  and  corresponds  to  the  time  of  beam  termination.  The  NO 
vibrational  bands  are  greatly  obscured  by  the  sharp  spiked  features,  which  . 
have  widths  equal  to  the  experimental  re  solution.  Since  unre  solved  vibration- 
rotational  bands  would  have  spectral  widths  much  greater  than  this, the  features 
were  suspected  to  be  either  due  to  atomic  lines  or  individual  rotational  lines 
of  a  molecular  vibrational  band. 

Under  certain  conditions  the  unidentified  spectral  features  were  found 
to  persist  long  after  NO  fluorescence  had  decayed.  Figure  13  shows  a  spectrum 
taken  in  the  same  experiment  as  Fig.  12,  but  at  a  time  corresponding  to 
3.4  msec  after  beam  termination.  Some  residual  NO  radiation  from  the 
slowest  decaying  observable  level  (v=2)  might  be  present  between  3650  and 
3770  cm"*.  Most  striking  .however,  is  the  appearance  of  the  unidentified 
spectral  features.  Indeed  some  weaker  features  attributed  to  noise  in  the 
other  runs  are  observed  to  be  reproducible  with  time  under  these  conditions. 
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Note  that  the  relative  intensity  of  the  lines  appears  unchanged  from  early 
time  observations,  see  comparison  of  Figs.  12  and  13.  These  long  time 
observations  provided  a  much  clearer  picture  of  the  unidentified  radiator. 
Specifically  the  strong  features  at  ~  3560  cm  ^  was  noted  to  agree  well  with 
the  frequency  of  the  Q  branch  of  the  OH(v=l-»0)  transition.  The  free  radical 
OH  was  considered  a  likely  radiation  source  inasmuch  as  water  was  intro¬ 
duced  into  the  system  with  the  room  air  and  could  have  a  long  residence 
time  on  the  chamber  surfaces  for  typical  system  pumping  rates. 

In  order  to  verify  the  tentative  assignment,  a  synthetic  spectrum  of 

,  2 

the  OH  molecule  was  generated.  The  OH  free  radical  has  an  inverted  II 
ground  electronic  state  with  strong  coupling  between  spin  and  angular  momentum 
giving  rise  to  two  total  angular  momentum  manifolds  of  states  Q=l/2  and  3/2.  - 
Transitions  between  similar  states  within  these  manifolds  are  not  degenerate, 
and  should  be  easily  observable.  Coupling  of  rotation  to  the  total  angular 
momentum  is  a  complex  function  of  the  rotational  and  vibrational  energy  and 
a  vibrating  rotor- anharmonic  oscillator  model  of  the  spectrum  is  not  suf¬ 
ficient.  Detailed  predictions  of  the  OH  spectral  line  positions  and  strengths 

14 

have  been  performed  by  Mies  and  shown  to  be  in  good  agreement  with 
experimental  spectra.  Mies*  calculated  frequencies  and  Einstein  coeffi¬ 
cients  were  used  to  generate  a  synthetic  OH'  spectrum  using  a  sine  slit  function 
and  assuming  a  300°K  rotaticnal  temperature.  A  comparison  between  this 
synthesized  spectrum  and  the  long  time  spectrum  of  Fig.  13  is  presented  in 
Fig.  14.  (Note  that  the  synthetic  spectrum  exhibits  negative  values  of 
intensity.  This  is  a  result  of  convoluting  the  predicted  OH  spectrum  over  a 
sine  slit  function  which  displays  prominent  negative  sidelobes.  This  effect 
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Figure  14:  Comparison  of  synthetic  spectrum  of  OH  with  spectral  features 
of  Fig.  13.  All  significant  features  (except  3270  cm'  line)  are 
matched  by  OH  predictions. 
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is  an  artifact  of  the  Fourier  transform.  )  It  can  be  seen  that  the  prominent 
spectral  features  match  OH  vibration/rotation  lines.  In  addition,  excellent 
agreement  occurs  between  the  spectral  positions  of  nearly  all  of  the  weaker 
features  and  OH  lines  arising  from  high  J  levels  and  the  weaker,  fl 

1  /  2 

transitions.  Note  that  in  the  synthetic  spectrum  the  relative  populations  of  the 
two  spin  states  and  the  rotational  levels  are  assumed  to  be  in  the  thermal 
equilibrium  with  the  gas  kinetic  temperature.  The  observed  relative  OH 
vibrational  populations  can  be  readily  determined  from  the  measured  Q 
branch  intensities  and  it  has  found  that  the  relative  populations  of  OH  v=2 
and  v=l  are  nearly  equal  at  beam  termination  (n,/n  =0.  8)  for  all  the  runs 
in  which  they  could  be  observed  (including  trace  oxygen,  and  fast  and  slow 
flow  cylinder  gas  and  room  air  runs.) 

Upon  backchecking  it  was  found  that  the  OH  spectrum  was  not  observed 
in  early  runs  made  using  cylinder  gases  prior  to  that  admission  of  room  air 
into  the  chamber.  Thus  it  appears  that  once  water  vapor  was  -omitted  it 
strongly  adhered  to  the  chamber  surfaces  and  ultimately  gave  rise  to  the  OH 
fundamental  band  radiation  even  in  "dry"  air  mixtures.  Additionally,  not  all 
the  ambient  water  vapor  is  removed  by  the  trap  in  the  filtered  room  air 
measurements. 

Due  to  the  slow  time  response  of  the  PbS  detector,  only  estimates  of 

the  quenching  rates  of  OH  by  O 2  and  could  be  made.  The  total  OH  quenching 

15 

rate  constant  by  an  air  mixture  of  N2/02  has  been  measured  as  1.  5  x 
-14  3 

10  cm  /sec,  dominated  by  O 2  quenching.  Observations  from  the  present  * 
work  support  O 2  quenching  of  OH  as  faster  than  N2»  but  slower  than  O^ 
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quenching  of  NO.  Quenching  of  OH  by  N^  appears  to  be  significantly  faster 

than  quenching  of  NO  by  N^,  and  is  estimated  to  have  a  rate  constant  of  order 
-15  3 

1.  5  x  10  cm  / sec  which  is  approximately  10%  of  the  rate  constant  for 
quenching  by  O These  observations  arise  from  the  analysis  of  the  low 
pressure  trace  oxygen  runs.  The  slow  decay  of  the  OH  radiation  in  Figs.  12 
and  13  (P=100  torr)  can  possibly  be  explained  by  vibrational  pumping  from 
N^vaZ)  which  would  not  make  as  large  a  contribution  in  the  low  pressure 
runs  (a  more  detailed  investigation  is  obviously  required).  Observations 
of  the  decay  OH(v=2)  are  limited  and  estimates  of  the  quenching  rates  are  not 
possible.  It  should  be  emphasized  that  the  present  study  was  not  designed 
to  study  OH  kinetics  and  the  above  observations  must  be  considerd  qualita¬ 
tive  rather  than  quantitative.  For  example  the  possible  chemical  destruction 
of  OH,  through  reactions  with  beam  produced  species,  has  not  been  considered. 
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A  number  of  interesting  side  features  which  could  be  studied  in  more 
detail  have  been  observed  in  this  data.  These  include:  (a)  the  observation 
of  OH  fluorescence  and  decay  which  might  be  utilized  to  ascertain  fundamental 
quenching  rate  constants  for  vib rationally  excited  OH,  (b)  NO  fluorescence 
data  sufficiently  resolved  to  allow  determination  of  vibrational  level  depen¬ 
dent  quenching  rate  constants  and  relative  fundamental  to  first  overtone 
band  Einstein  coefficients,  and  (c)  the  possibility  of  using  trace  quantities 
of  CO?  to  determine  relative  populations  of  N-O  and  vibrationally  excited  N-. 
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